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Abstract 


Multivariate morphometric analyses of cranial and 
post-cranial skeletal data from historical, Holo- 
cene, late Pleistocene and living remnant popu- 
lations of bison were used to study geographic 
variation and subspeciation of Recent bison. Par- 
ticular consideration was given to Bison bison 
athabascae and the taxonomic status of the rem- 
nant populations of presumed pure wood bison 
from the Nyarling herd. External characters 
(standard measurements, weight, outline shape 
and pelage characters) were compared in living 
forms. 

Geographic variation in historic populations of 
the North American bison appears to have been 
clinal in nature over most of the distributional 
range in a predominantly south — north direction. 
However, a morphological discontinuity existed 
in the northwestern part of the range, coinciding 
approximately with the grassland — boreal for- 
est ecotone. This phenotypic discontinuity between 
grassland and woodland populations fully justi- 
fies recognition of the two current subspecies B. b. 
bison and B. b. athabascae. The historical distri- 
bution of the two subspecies was reconstructed 
using information from specimens, historical 
reports and toponyms. 


vl 


Comparisons of modern American and Euro- 
pean bison with late Pleistocene — Holocene bison 
Suggest, the existence of one variable chrono- 
species in the Holarctic region. The presence of 
bison closely similar to B. b. occidentalis and B. 
b. athabascae in eastern Siberia during the late 
Pleistocene-early Holocene demonstrates the exis- 
tence of a Beringian population ancestral to living 
bison. It is probable that modern North American 
bison originated from the Beringian form and a 
southern form (B. b. antiquus). Subsequent evo- 
lution, coinciding with the expansion of the grass- 
lands gave rise to B. b. bison. B. b. bison is the 
more derived of the current subspecies, whereas 
B. b. athabascae retains more of the ancestral 
characters, resembling B. b. occidentalis to a 
greater degree. Comparisons of cranial and, post- 
cranial elements, as well as those of external 
characters show that the Nyarling River bison are 
close to the original B. b. athabascae. However, 
the evidence also suggests that some gene flow 
between this herd and hybrid herds elsewhere in 
Wood Buffalo National Park did occur. 


Résumé 


Des analyses morphomeétriques multivariables de 
données craniennes et squelettiques provenant 
d’échantillons de spécimens de bisons datant de 
l’Holocéne et du Pléistocéne, ainsi que de popula- 
tions reliques ont été entreprises afin de clarifier 
des problémes relatifs a la variation géographique 
et a la sous-spéciation des bisons contemporains. 

On a considéré en particulier le Bison bison 
athabascae et la position taxinomique des bisons 
des foréts présumés pures qui ont survécu dans la 
région du fleuve Nyarling. Des caractéres exter- 
nes (mensurations types, poids, contour du corps 
et pelage) ont été comparés chez les formes 
vivantes. 

La variation géographique chez les populations 
historiques de |’ Amérique du Nord était apparem- 
ment clinale dans la plus grande partie de I|’aire 
de répartition, surtout dans une direction nord- 
sud. Cependant une discontinuité morphologique 
existait dans la région du nord-ouest, adjacente 
a l’écotone entre la steppe et la forét boréale. Cette 
discontinuité phénotypique entre les populations 
des steppes et de la forét justifie les sous-espeéces 
actuelles B. b. bison et B. b. athabascae. La répar- 
tition historique des deux sous-espéces a été 
reconstruite en utilisant les endroits ot les spéci- 
mens ont été trouvés, des mentions historiques et 
des toponymes. 


vii 


Des comparaisons entre les bisons actuels 
d’Amérique et de |’ Europe et des bisons de la fin 
du Pléistocéne ou du début de l’Holocéne suge- 
rent qu’il y avait une espéce variable dans la région 
holarctique. L’existence d’un bison proche de 
B. b. occidentalis et de B. b. athabascae en Sibérie 
orientale vers la fin du Pléistocéne ou le début de 
l’Holocéne démontre qu’il y avait une population 
ancestrale du bison actuel en Béringie. 

Le bison actuel d’Amérique du Nord a proba- 
blement eu comme ancétre cette forme béringienne 
ainsi qu’une forme méridionale (B. b. antiquus). 
L’évolution durant l’expansion des steppes, a plus 
tard donné naissance au B. b. bison. Ce dernier 
est le plus évolué des sous-espéces actuelles, tan- 
dis que le B. b. athabascae a conservé plus de 
caractéristiques ancestrales plus marquées et res- 
semble plus au B. b. occidentalis. Une comparai- 
son des caractéristiques squelettiques et externes 
démontre que le bison du fleuve Nyarling est pro- 
che de B. b. athabascae originaire de cette région. 
Les données indiquent en outre qu’il y avait un 
échange de genes entre ce troupeau et les autres 
troupeaux hybrides, ailleurs dans le parc national 
de Wood Buffalo. 
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Introduction 


More than half a century ago, 6673 plains bison 
(Bison bison bison) were introduced to Wood 
Buffalo National Park, the range of the last sur- 
viving 1500 to 2000 wood bison (B. b. athabascae) 
(Seibert, 1925). Since then, it has been generally 
assumed that the wood bison had disappeared as 
a distinct subspecies as a result of interbreeding 
with the more numerous plains bison. The pheno- 
type of the plains bison is now said to predomi- 
nate in the Park (Fuller 1962). In 1957, however, 
an isolated herd of approximately 200 bison was 
discovered in the northwestern part of Wood 
Buffalo National Park, which was assumed to be 
a remnant of pure wood bison. This assumption 
was based on the observation that the intervening 
terrain, consisting of impassable and unsuitable 
habitat, formed a barrier to the movement of 
bison. Subsequent comparison of three specimens 
from this population with specimens of B. b. 
athabascae \ed to the conclusion that they were 
indeed wood bison (Banfield and Novakowski 
1960). Shortly after the discovery of the remnant 
herd, the Canadian Wildlife Service undertook a 
program to save the herd from possible genetic 
inundation and.to secure the survival of the 
subspecies (Reynolds 1979). 

In the context of the Wood Bison Rehabilita- 
tion Program, the more extensive, in-depth anal- 
ysis of the taxonomic status of the remnant herds 
of presumed pure wood bison was required, and 
the present study was undertaken to meet that 
need. 

Although the emphasis in the present study is 
thus on B. b. athabascae and the identity and tax- 
onomic status of the remnant herds of presumed 
wood bison, the analysis required the examination 
of the problem in the broader perspective of geo- 
graphic variation and the systematics of historical 
bison populations and, to some extent, earlier 
forms of bison as well. 

The near extermination of the bison in North 
America in the latter part of the nineteenth century 
precluded a comprehensive analysis of geographic 
variation in modern bison on this continent. The 
two North American subspecies recognized at 
present may therefore, not adequately reflect the 
geographic variation that once existed. It appears 
indeed likely that this species with its vast historical 
distribution may have had more distinct geo- 
graphic forms than are now recognized. Some 
authors (Heck 1936, Mohr 1952, Krumbiegel 1980) 


have suggested, chiefly on the basis of historical 
illustrations and photographs, that the plains 
bison comprised at least two distinct geographic 
forms. 

The recognition of the existence of a distinct 
form of bison in northern Canada can be traced 
back to the latter part of the last century (Allen 
1876, Seton 1886) and culminated in the formal 
description of the subspecies B. b. athabascae by 
Rhoads. Rhoads’ (1897) description was based on 
a single specimen which was not examined directly 
by him. The lack of a clear and unambiguous con- 
cept of the subspecies, resulting from the lack of 
critical comparisons, caused widely varying opin- 
ions on the taxonomic status of B. b. athabascae 
to persist to our day. Soper (1941: p.375) summed 
up the situation succinctly when he wrote: ‘‘Some 
regard the differences [between wood and plains 
bison] as but slight and of little significance at 
most a superficial divergence springing from 
environment; others again support Rhoads in this 
definite separation; while a third group leans to 
recognition of its status as a distinct species.”’ 

While a majority of mammalogists has sub- 
scribed to the view that athabascae is at least 
subspecifically distinct (e.g. Allen 1876, Seton 
1886, Rhoads 1897, Rowan 1929, Skinner and 
Kaisen 1947, McDonald 1981), there have always 
been dissenters. For example, Hornaday (1889), 
Graham (1923) and Garretson (1938) would fit in 
Soper’s first group, while Flerov (1979) represents 
a recent proponent of the third. 

The lack of first hand knowledge of the living 
animal-only Seton and Rowan among zoologists 
could claim to have seen the plains and wood bison 
alive — and the small number of specimens avail- 
able for study were undoubtedly at the root of this 
diversity of opinions. Quantitative analyses of 
B. b. athabascae museum material, primarily 
crania, were not undertaken until 1947 by Skinner 
and Kaisen and 1978 by McDonald. The ill- 
conceived introduction of the plains bison into the 
last refuge of the wood bison apparently doomed 
the latter’s survival as a distinct subspecies before 
a detailed study of all aspects of this bison could 
have begun. The historic discovery in 1957 raised 
the hope that not all was lost. The aim of the 
present study is to determine to what extent the 
remnant saved represents the original wood bison 
that once roamed the area of what is now Wood 
Buffalo National Park. 


The principal questions addressed in the present 
study are: 

1. What was the extent and nature of geographic 
variation in North American Recent bison? 
2. What are the morphometric relationships of 

historical North American forms to earlier 

Holocene and the late Pleistocene forms, and 

what are the implications for the origin of 

Recent bison? 

3. How are recent New and Old World bison 
related? 

4. What is the taxonomic status of present remnant 
herds in Canada? 

The first three of these questions have been 
addressed before. This study is, therefore, chiefly 
a reassessment, using a somewhat different 
approach and different analyses. An analysis of 
the evolution of bison since the early Pleistocene 
lies outside the scope of the present study. For a 
current treatment of this subject see Weston, 
Shackleton and Harington’s (in preparation) com- 
prehensive analysis of fossil and recent material 
from Eurasia and North America. In the context 
of the present study, it suffices to glance back to 
the earlier Holocene and late Pleistocene in order 
to gain some insight into the possible origin of 
contemporary forms of bison. 

The fourth question focuses on the survival of 
the wood bison as a distinct subspecies. Aside 
from its scientific interest, the answer to this ques- 
tion may have practical implications for bison 
conservation and management. 

The effects of interbreeding of the plains and 
wood bison on the phenotypic variation of the 
present herds in Wood Buffalo National Park 
could, unfortunately, not be fully assessed in the 
present study because of a lack of specimens. 
Further work on this population will be under- 
taken. 


Materials and Methods 


The analyses in this study are based on data 
derived from 528 skulls; 143 mandibles; 36 post- 
cranial skeletons; a smaller number of skins as well 
as on external measurements and weights of fresh 
carcases and information on external characters 
from live animals gathered in the course of this 
study. 

Pertinent information on the specimens used 
and number of specimens in each sample are given 


in Appendix 1. 

The specimens were divided a priori for pur- 
poses of analysis into a number of subsamples, 
on the basis of sex, geographic origin and era as 
follows: 


Subsamples Representing Living Remnant 
Populations. 


Subsample 1. Descendants of the Allard-Pablo 
herd (AP). This sample consists chiefly of speci- 
mens from Elk Island National Park, with smaller 
numbers of specimens from Wainwright (former 
Buffalo National Park) and Riding Mountain 
National Park. 

Subsample 2. Wood Buffalo National Park 
(WBNP), consisting of specimens collected 
20 years or more after the final release of AP stock 
and presumed to represent the mixed descendants 
of it and athabascae. Only an inadequate sample 
of mainly immature animals could be gathered, 
precluding a complete analysis of this population. 
Subsample 3. Snyder Expedition (SE). This small 
sample consists of specimens collected in Wood 
Buffalo National Park in 1934-35 (Goodwin 
1935). The aim of the expedition was to collect 
specimens of wood bison for the American 
Museum of Natural History and the National 
Museum of Canada. Selection was on the basis of 
external appearance by implication, but as these 
animals were collected 6 or 7 years after the last 
introduction of plains bison into the park and 
none of the collectors had prior first-hand knowl- 
edge of pure B. b. athabascae, the taxonomic iden- 
tity of these specimens is in doubt. Because of the 
definite bias used in selecting the specimens, this 
sample cannot be used to represent the hybrid 
population in Wood Buffalo National Park as 
McDonald (1978, 1981) has done. 

Subsample 4. Nyarling River herd (NR), a sample 
comprising specimens of presumed pure B. b. 
athabascae discovered in 1957 and their descend- 
ants. Small numbers of animals transferred from 
this herd to other areas (18 animals to the Mack- 
enzie Bison Sanctuary in 1963 and 23 animals to 
Elk Island National Park in 1965) (Reynolds 
1979), formed the foundation of all present wood 
bison herds. 

Subsample 5. bonasus, a sample representing the 
European bison (B. b. bonasus and a few B. Db. 
caucasicus) and consisting chiefly of craniometric 
data made available by Dr. H. Bohlken. 


Subsamples Representing Historical Populations 


Subsample 6. athabascae, consisting of specimens 
of B. b. athabascae collected prior to the 1925-28 
introductions of plains bison to Wood Buffalo 
National Park, or those collected during or shortly 
after that period, that could be unequivocally 
assigned to B. b. athabascae on the basis of age. 
The buik of the 6673 plains bison introduced con- 
sisted of yearlings (72.3%) and two-year-olds 
(22.7%) the remainder of three-year-olds (Soper 
1941). Only one true B. b. athabascae female 
specimen (NMC 10405) was available for study. 
McDonald’s (1978, 1981) rejection of this speci- 
men as a true B. b. athabascae appears to be with- 
out foundation. The female’s age was determined 
on the basis of the number of incremental lines 
in her tooth cementum as 10+ years. As she was 
collected in March 1928, she was therefore at least 
7 years old when the first introductions were made 
in 1925. 

Subsample 7. Canadian Prairies (CP) consists of 
specimens representing the original historical 
population of bison from the Canadian prairie 
provinces. 

Subsample 8. Canadian Rocky Mountains (CRM), 
a sample consisting of historical specimens, sur- 
face and near-surface finds, from Banff, Jasper 
and Waterton Lakes National Parks. 
Subsample 9. Peace River region (PR) consisting 
of skulls and partial skulls of historical specimens 
from the Aspen-parklands of the Peace River 
region of Alberta and British Columbia. This area 
is of particular interest with regard to the histori- 
cal distribution and geographic variation of bison. 
As there were few specimens from this area in 
collections, a publicity campaign was mounted 
through local newspapers to solicit information 
about specimens picked up locally and in private 
ownership. In this manner additional specimens 
from the area were located and examined. 
Subsample 10. Central Great Plains (CGP) con- 
sisting of specimens representing historical popu- 
lations from the central Great Plains in the United 
States, including specimens from Iowa, Kansas, 
Montana, Nebraska, South Dakota and eastern 
Wyoming. 

Subsample 11. United States Rocky Mountain 
(USRM), consisting of historical specimens from 
the Rocky Mountains and areas to the west, 
including specimens from Colorado, Montana, 
western Wyoming, and Oregon. 


Subsample 12. Southern Great Plains (SGP) con- 
sisting of specimens representing the historical 
population in Texas. 

Subsample 13. Eastern United States (EUS) con- 
sisting of a small sample from Kentucky. 


Subsamples Representing Earlier Holocene and 
Late Pleistocene Populations 


Subsample 14. occidentalis representing early to 
mid-Holocene North American bison. 
Subsample 15. antiquus representing late Pleisto- 
cene to Holocene North American bison. 
Subsample 16. priscus representing late Pleisto- 
cene Beringian bison from Alaska. 


In addition to the above samples data from 
individual skulls of Holocene or late Pleistocene 
bison were used (see Appendix 1). 

The subsamples are referred to in the text by 
the taxonomic designations or the acronyms listed 
above. Cranial measurements were in some 
instances partially (CGP; USRM: occidentalis), or 
wholly (SGP; EUS; antiquus; priscus) derived 
from McDonald (1978). For further information 
On specimens used see Appendix 1. 

The cranial and mandibular measurements used 
are those of Skinner and Kaisen (1947), and 
Shackleton et al. (1975) with a few additional 
ones (see Figure 1). Post-cranial measurements 
and external measurements taken are shown in 
Figures 2 and 3 respectively. The abbreviations of 
the measurements used are listed in figures 1 to 3. 

The numerical data were analyzed using univar- 
late and multivariate methods. Basic statistics were 
calculated for all measurements and the sig- 
nificance of differences between sample means 
tested with Student’s ¢ test. Use of single varia- 
bles for discriminating between subspecies or sexes 
was evaluated using the method of Lubischew 
(1962). Cranial, post-cranial and external varia- 
bles were also subjected to principal component 
analysis (PCA) (in cases involving individual speci- 
mens) or canonical variates analysis (CVA) (of a 
priori group) and distance analysis (Blackith and 
Reyment 1971, Neff and Marcus 1980, Davies 
1971). Conformity of the data to the assumptions 
of multivariate normality and equal variance 
covariance structure were not tested. According 
to Sokal and Sneath (1973: p. 127) ‘‘. . .consider- 
able robustness to violation of these assumptions 


Figure 1. Cranial and mandibular measurements used. Numbers correspond to those of Skinner and Kaisen (1947) and Shackleton 
et al. (1975). Cranial measurements: 1) Spread of horn cores, tip-to-tip (SHTT); 3) Core length on upper curve, tip to burr (CLUC); 
4) Core length on lower curve, tip to burr (CLLC); 5) Length, tip of core to upper base of burr (TB); 6) Vertical diameter of 
horn core at right angles to longitudinal axis (VD); 7) Circumference of horn-core at right angle to longitudinal axis (CHC); 8 Greatest 
width at auditory openings (GWA); 9) Width of condyles (CW); 12) transverse diameter of core at right angle to longitudinal 
axis (TD); 14) Width of cranium between horn-cores and orbits (WHCO); 15) Greatest postorbital width (GPW); 16) Anterior 
orbital width at notch (AOW); 17) Width of skull at masseteric processes above M1 (WMP); 18 a and b) Rostral width at maxillary- 
premaxillary suture (RW1, RW2); 19) P2-M3 alveolar length (P2M3L); 20) M1-M3 alveolar length (M1M3L); 23) Width of the 
occipital constriction (OC); FP basilar length; ON length occipital crest to nasofrontal suture. Mandibular measurements: 
7) ml-m3, alveolar length (m1m3L); 9) height of the coronoid process (HCP); 10) Width of coronoid process (WCP); 11) Width 
of articulating process at mandibular notch (WAP); 12) Posterior of m3 to mandibular foramen (Pm3F); 13) Posterior of m3 
to posterior edge of articulating process (Pm3AP); 14) Width of ramus at posterior of ml (WR); 16) Width of mandibular condyle 
(WMC); 17) m2 alveolar to base height (mandibular depth) (MD); 21) pm2-m1 (premolar length) (PL); 22) substitute total length 
(STL); 24) diastema width (DW). 
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Figure 2. Measurements of the postcranial skeleton. 1) Length of scapula (SCAPL); 2) Width of scapula (SCAPW); 3) Length 
of humerus (HUML); 4) Least width of humerus (HUMW); 5) Length of radius (and ulna) (RADL); 6) Proximal width of radius 
(RADW1); 7) Distal width of radius (RADW2); 8) Length of metacarpal (MCL); 9) Proximal width of metacarpal (MCW1); 10) 
Distal width of metacarpal (MCW2); 11) Length of femur (FEML); 12) Proximal width of femur (FEMW1); 13) Distal width 
of femur (FEMW2); 14) Length of tibia (TIBL); 15) Proximal width of tibia (TIBW1); 16) Distal width of tibia (TIBW2); 17) Length 
of metatarsal (MTL); 18) Proximal width of metatarsal (MTW1); 19) Distal width of metatarsal (MTW2); 20) Length of pelvis 
(PELV); 21) Width of ilium (ILW); 22) Width at acetabulum (ACETW); 23) Width at ischium (ISCHW); 24) Length of neural 
spine of vertebrae, last cervical, thoracic, and lumbar (C7, T1-14, L1-5). 
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Figure 3. A. External measurements: 1) Total length (TL); 2) Length of the tail a) to last caudal vertebra (T); b) to end of tassel. 
3) Height of the shoulder (SHHT); 4) Rump height (RHT); 5) Length of hindfoot (HF). B. Terminology used for pelage characters, 
the broken line outlines the area of the cape. C. Standardization for outline-shape analysis. The divisions on the horizontal axis 
are in centimeters. 


has been demonstrated’’. Moreover the descrip- 

tive use of multivariate techniques ‘‘. . . as a data- 

analytic procedure. . . to ordinate and efficiently 
describe multivariate observations’’ (Albrecht 

1980: p. 681) is valid. Discriminant function 

analysis (DFA) for two groups for use in assigning 

unknowns to B. b. bison or B. b. athabascae was 
also used (Appendix 2). Calculations were per- 
formed on a Hewlett-Packard 9831 A desk-top 

computer, using programs from Davies (1971) 

adapted by D. McAllister. 

During preliminary analysis the initial number 
of cranial variables, 34 in all, was reduced to a 
smaller number (7 to 11) in subsequent analyses. 
Deletion of variables was based on redundancy 
reflected by high correlations between variables 
and poor representation of a variable in the B. b. 
athabascae sample in order to maximize sample 
size for this taxon. The initial number of mandibu- 
lar variables (23) was likewise reduced to 15 in the 
final analyses. 

In the craniometric analyses data from adult 
animals only were used. Assignment of skulls to 
age classes was based on the degree of fusion of 
the parietofrontal and frontal sutures in males 
(Shackleton et al. 1975) and toothwear in both 
sexes (Skinner and Kaisen 1947; Fuller 1959b; 
Frison and Reher 1970). 

To elucidate the questions posed in the introduc- 
tion the following analyses were undertaken: 

1. The morphometric relationships of historical 
populations (subsamples 6-13) were assessed 
using CVA and distance analysis of cranio- 
metric data. CVA obtains a low dimensional 
display of the data and maximizes between 
group variation relative to within group varia- 
tion (Campbell and Atchley 1981). Plots of the 
canonical variates scores in two and three 
dimensions were produced to display the scat- 
ter of individuals and the distance between 
group centroids. To assess the contribution of 
the original variables to the canonical variates 
the vectors scaled to have a length of two stand- 
ard deviations were plotted on the graph 
(Jolicoeur 1959). The canonical variables were 
used to calculate a matrix of Mahalanobis dis- 
tances (Pythagorean distance in standardized 
canonical variate space). To facilitate interpre- 
tation of the relationships among groups in 
three-dimensional canonical variate space and 
to reveal possible distortions of distances, mini- 
mum spanning trees (MST) (Sneath and Sokal 
1973) were computed and projected on the 
graph. The distance (D) relationships among the 


samples were in some instances further summa- 
rized by a distance dendrogram derived by 
unweighted pair group cluster analysis 
(UPGMA) (Sneath and Sokal 1973) of the 
distance half-matrix. 

The 11 cranial variables used in the analysis 
included 6 horn measurements, which showed 
high intra-locality correlations suggesting com- 
mon epigenetic control and redundancy of 
information (Thorpe 1976). This set of variables 
was used, however, so that some of McDonald’s 
(1978) data could be included in the analysis. 

Different sets of variables, in which redun- 
dancy had been reduced by selecting those 
variables with low intra-locality correlations 
(Thorpe 1976), were used to analyze the mor- 
phometric relationships of northern populations 
CP, CGP, CRM, PR and B. b. athabascae fur- 
ther. To determine the extent of the historical 
distribution and area of possible intergradation 
of the two forms, all available specimens from 
critical areas, many being partial skulls, were 
used. Histograms of variables most frequently 
preserved in historical specimens were con- 
structed (Figure 4) to aid in selecting combina- 
tions of characters useful in assigning unknown 
skulls to the correct taxon with a degree of prob- 
ability. Discriminant function analysis, using 
combinations of these variables varying from 
11 to as few as 3 were used to assign unknown 
specimens from critical areas (Appendix 2). 
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Figure 4. Histogram showing the frequency of intact measur- 
able variables in old skulls found on or near the surface of 
localities in the historic range of the bison. 


2 .The morphometric relationships of Recent sub- 
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species to earlier Holocene and late Pleistocene 
forms (occidentalis, antiquus, priscus) was 
analyzed in the same manner as that of the 
historical populations and by PCA. 


. The relationship of modern New and Old 


World bison was assessed by CVA of 11 cranial 
variables. A somewhat different set of cranial 
variables, corresponding to that of Bohiken 
(1967) was used. 


. The taxonomic status of living remnant herds 


in Canada was determined using CVA of cra- 
niometric data and PCA of postcranial variables 
and external measurements and weights. In 
addition, outline shape and pelage chracteris- 
tics of NR bison were compared with those of 
plains bison. 

The PCA of post-cranial variables was done 
on the variance-covariance matrix and correla- 
tion matrix of the data. External measurements 
and weights were transformed into logarithms 
(logjo) prior to PCA of the variance-covariance 
matrix of the data. 

Differences in outline shape and pelage 
characters of American plains bison and Euro- 
pean bison are well known and discussed in 
detail by Mohr (1939, 1952) and Flerov (1932), 
and summarized by Haltenorth (1963). Infor- 
mation on the external characteristics of the 
original B. b. athabascae, on the other hand, 
is poorly documented. Aside from a few 
museum skins, which give information on pel- 
age characters, there is very little objective infor- 
mation on the external characteristics of B. b. 
athabascae. Early reports of eye witnesses 
appear to agree on two characters in which B. Db. 
athabascae differed from B. b. bison, i.e. size 
(larger) and colour (darker). Existing photo- 
graphs (e.g. Graham 1923, Kitto 1924) predat- 
ing the 1925-28 introductions are generally 
inadequate for a detailed comparison of the two 
forms. In the absence of a B. b. athabascae 
reference sample for external characters, NR 
specimens were compared directly to AP speci- 
mens and to a lesser extent also to plains bison 
from different herds. 

Outline tracings were obtained from photo- 
graphs of bison in profile. Outlines were stand- 
ardized so that each had an arbitrary length of 
16 cm along a horizontal axis from the base of 
the tail to the plane of the nose (Figure 3). The 
two orthogonal axes superimposed on the out- 
line allow for a quantitative comparison of 


differences in outline shape. Perpendicular dis- 
tances of contour points (at intervals of 1 cm 
on the X axis) to these axes could thus be 
recorded and compared. 

Pelage characters in NR and AP bison were 
compared using visual character analysis, com- 
parison of frequencies of occurrence of certain 
characters and a character index (Mayr 1969). 
A number of characters which were found to 
differ in NR and AP samples after preliminary 
examination were subdivided and ranked into 
two or more character states, whose frequency 
of occurrence in each group could then be deter- 
mined and compared. Characters showing the 
least amount of overlap were then assigned a 
numerical value, and could thus be used in a 
character index. The typical or prevailing 
character state of AP bison was designated 
as 1. The character state encountered in the NR 
sample most unlike that found in the AP sam- 
ple was designated 0. Where intermediate 
character states are recognized they were 
assigned intermediate values (see Appendix 3). 
In interpreting the data it was assumed that 
character states most unlike those found in AP 
bison or absent in AP bison were characteristic 
of B. b. athabascae. Where possible, confirma- 
tion was sought from independent sources. With 
the addition of one outline shape character the 
character index was used to score individual 
animals to obtain an index of their position on 
an axis representing external phenotypic varia- 
tion from B. b. bison (100%) to B. b. atha- 
bascae (0%). 


Results 


1. CVA of Craniometric Data from Historical 
Populations of North American Bison. 


The results of the CVA and distance analysis of 
samples representing historical populations from 
different parts of the former distributional range 
based on 11 cranial variables are summarized in 
Figures 5, 6 and 7. 

Examination of the three-dimensional graph 
(Figure 5) clearly shows the isolated position of 


the athabascae centroid relative to the other forms. 


The centroids for the populations, excluding 
athabascae which may be collectively referred to 
as plains bison, are separated by relatively small 
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Figure 5. Projection of the centroids for eight samples of historical North American bison populations onto the first three canonical 
axes, based on six horn measurements and five cranial measurements. The centroids are connected by a minimum spanning tree. 
1) Eastern U.S. 2) Southern Great Plains; 3) U.S. Rocky Mountains; 4) Canadian Rocky Mountains; 5) Peace River area: 6) Canadian 
prairies; 7) Central Great Plains; 8) B. b. athabascae. 
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Figure 6. Projection of individual specimens and centroids of eight samples onto the first two canonical axes. The projection 
of character vectors shows the relative contribution of the original variables to the canonical variates. The vector lengths equal 
two standard deviations. The ellipses represent the 95% confidence region of the centroids. 1 ON; 2 VD; 3 TD; 4 GWA; 5 CW; 
wewtiCcO: 7 GPW; 8 SHTT; 9 CLUC; 10 TB; 11 CHC. 
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Figure 7. Dendrogram summarizing the generalized distances 
between historical populations of North American bison in 
canonical variate space (cf. Figure 5.) The cophenetic 
correlation coefficient is 0.87. 
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distances (D) from their nearest neighbors in 
canonical variate space as the MST shows. The 
distances linking the plains bison centroids vary 
from 0.76 (between CP and PR) to 2.78 (between 
EUS and SGP), whereas athabascae and its nearest 
neighbor PR are considerably more distant from 
each other (D = 4.25). In addition to their prox- 
imity to one another the centroids of the plains 
bison samples are distributed along a diagonal axis 
running from the left back upper corner down 
toward the center of the three-dimensional space. 
The centroid of athabascae, on the other hand, 
is situated well to the right of this axis. This axis 
represents clinal variation in the plains bison in 
an approximately north-south direction. 

The projection of the canonical variates scores 
for individual specimens and their group centroids 
on the first and second canonical axes shows that 
the plains bison samples form a series of overlap- 
ping clusters (Figure 6). The clinal variation within 
the plains bison, referred to above, appears to be 
predominantly that of overall size (Table 1), with 
the smaller forms in the upper left corner and the 
larger ones lower down towards the center. 

The relative contribution of the original varia- 
bles to the canonical variates is shown by the 
vector diagram on the plot of canonical variates 
(Figure 6). The vector separating the plains bison 
centroids from each other is GWA (no. 4), a meas- 
urement significantly correlated (r=0.81 P<.02) 
with overall size (expressed as the mean of means 
of all measurements for each group). The vector 
for GWA runs parallel to the axis of clinal varia- 
tion in the plains bison suggesting that the 
variation is essentially one of size. 

The athabascae cluster, on the other hand, does 
not overlap any of the plains bison clusters and 
is clearly divergent from the axis of clinal variation 
of the plains bison. The position of athabascae 
suggests a difference in shape of the athabascae 
crania compared to those of plains bison. The 
character vectors (Figure 6) indicate that athabas- 
cae differs from the plains bison sample in having 
greater TD, CLUC, GPW, WHCO. The distance 
relationship of the historical populations is shown 
in Figure 7. The distance dendrogram clearly 
shows the divergent position of athabascae. 

The morphometric relationship of the northern 
plains bison (CP, CGP, USRM, CRM and PR) 
and athabascae was subjected to further analysis 
using different sets of variables, ranging from 11 
to as few as 3, in order to elucidate variation in 
the grassland-boreal forest ecotone, the possible 
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Figure 8. Projection onto the first two canonical axes of individual specimens and centroids (connected by a minimum spanning 
tree) of six samples representing populations of northern plains bison and wood bison. CVA is based on 9 cranial and 2 horn 
characters. Lengths of character vectors equal two standard deviations. 1 ON; 2 VD; 3 TD; 4 GWA; 5 CW; 6 WHCO; 7 GPW 


8 P2M3L; 9 MIM3L; 10 OC; 11 WMP. 


zone of intergradation between wood and plains 
bison. 

The graph in Figure 8 summarizes the results 
of the CVA of 11 cranial variables, of which only 
two are horn measurements. The results are 
remarkably similar to the previous CVA of 11 
variables in which horncore measurements 
predominated. There is a clear separation between 
northern plains bison samples and athabascae. 
One specimen from the Peace River parklands and 
two specimens from the Canadian Rocky Moun- 
tains were scored as unknowns and plotted to 
determine their affinities. All three specimens fell 
within the scatter of the plains bison samples 
(Figure 8). The character vectors show that 
athabascae differs from the plains bison samples 
mainly in greater GWA, TD and smaller WMP 
and CW. 

A CVA using seven cranial variables with an 
increased number of specimens from the Peace 
River parklands produced a similar picture 
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(Figure 9), but with smaller inter-centroid dis- 
tances. The interesting result to note is the some- 
what intermediate position of PR in relation to 
athabascae and the plains bison samples. Similar 
patterns were observed in CVA’s using progres- 
sively fewer variables, with the main difference 
being the concomitant decrease of intercentroid 
distances and increasing overlap of clusters 
representing the geographic samples. As few as 
three variables still gave only approximately 4% 
overlap between athabascae and plains bison refer- 
ence samples and could therefore be used in 
assigning unknown, geographically intermediate 
specimens — a large number of them with only 
a few measurable variables (Figure 4), to the 
nearest reference sample. Partial specimens from 
the area between the main distributional ranges 
of plains and wood bison could thus be identified 
and used in delineating the contact zone between 
wood and plains bison. The geographic and mor- 
phometric relationships based on three variables 
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Figure 9. Morphometric relationship of northern plains bison and wood bison on the basis of a CVA of five cranial and two 
horn characters. 1 VD; 2 TD; 3 GWA; 4 CW; 5 WHCO; 6 GPW; 7 OC. 


(VD, TD, WHCO) of four geographic samples 
(CP, CRM, PR, athabascae) are shown in 
Figure 10. 

The use of three variables, although decreasing 
the discrimination, allowed most available speci- 
mens to be used in determining the area of con- 
tact or intergradation between plains and wood 
bison. The small number of measurements and 
decreased power of discrimination meant that 
large plains bison or small wood bison could be 
assigned incorrectly. For example a known large 
plains bison (NMC 8833) was scored as a wood 
bison based on the discriminant function using 
only three variables. However, a general trend and 
the affinities between geographic samples can be 
discerned from MST in Figure 10, with CRM and 
CP and PR and athabascae respectively being rela- 
tively close, while the two pairs are more distant 
from each other. It appears likely that the zone 
of contact was situated somewhere midway 
between CP and PR, judging from the relative 
proportions of bison and athabascae scores in the 
area. About 23% of the PR specimens scored on 
the bison side of the dividing point, compared to 
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about 50% of specimens from the boreal forest 
and aspen parklands section of Alberta separat- 
ing the Peace River parklands from the grasslands. 
Compared to the CVA using 7 variables and fewer 
specimens the most notable difference in the three 
variable analysis is the shift of the PR centroid to 
athabascae. The proximity of some PR specimens 
to CP suggests that the populations in the two 
areas were at least in occasional contact. 


2. CVA of Craniometric Data from Recent 
European and North American Bison. 


The morphometric relationship of crania of 
modern North American (historical samples) and 
European bison are shown in Figure 11. Prairie 
bison and European bison are separated from each 
other by the smallest distance. The main separa- 
tion between the two forms is along the second 
axis, with the most important contribution to dis- 
crimination coming from OC, AOW (greater in 
bonasus) and GWA, WMP, GPW (greater in 
bison). Discrimination between athabascae and the 
other two taxa is chiefly along the first axis with 
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Figure 10. Morphometric relationship of contiguous populations of bison from the Canadian prairies (CP), the Rocky Mountains 
(CRM), the Peace River area (PR) and from the area of what is now Wood Buffalo National Park (athabascae), based on three 
cranial characters. Solid dots represent specimens with values on the B. b. bison side of the mid point between the B. b. bison 
and B. b. athabascae centroid scores, hollow circles those with values on the athabascae side. The distances between groups show 
their mutual relationship and indicate that the most change occurs between CP and PR. Present vegetation zones are based on 


Rowe (1972). 


the greatest distance being that between atha- 
bascae and bonasus. Contribution to the discrimi- 
nation between athabascae and bonasus comes 
mainly from GWA, GPW, FP (greater in atha- 
bascae) and AOW, OC. Discrimination between 
athabascae and bison is contributed to mainly by 
GWA, GPW, FP and WMP (Figure 11). 
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3. Multivariate Craniometric Comparisons of 
Historical, Holocene and Late Pleistocene Bison. 


The results of the CVA of 11 cranial variables of 
historical, and early to mid-Holocene (occiden- 
talis, antiquus) and late Pleistocene forms from 
Alaska (priscus) are shown in Figure 12. The sepa- 


ration between smaller modern short-horned bison 
and larger long-horned Holocene and late Pleisto- 
cene bison is along the first axis. The contribution 
to the discrimination along this axis comes mainly 
from horn measurements (TD, TB, CLUC, 
SHTT) and to a lesser extent from cranial width 
measurements (WHCO), (cf. Table 1) The posi- 
tion of athabascae on the first axis is intermedi- 
ate with respect to the plains bison on the one hand 
and to occidentalis and antiquus on the other. 

The second axis shows most of the variation 
between the geographic samples of historical bison 
with the positions of athabascae and EUS on 
opposite extremes. The contribution to the dis- 
crimination along the second axis comes mainly 
from GWA, TD, TB. The high score of atha- 
bascae on this axis is owing in large part to the 
large dimension of GWA in this form. 

The distance relationships of the centroids are 
summarized by the MST in Figure 12 and the dis- 
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tance dendrogram in Figure 13. Examination of 
the latter clearly shows the remote position of 
priscus from the others, the intermediate position 
of athabascae between modern plains bison and 
occidentalis and antiquus, and, within the plains 
bison the relatively distant position of EUS. 

A well preserved skull of a bison from eastern 
Siberia (Moscow Paleontological Museum 835- 
624-39), believed to be of late Pleistocene or early 
Holocene age (Sher 1971), is of particular interest 
for the evolutionary history of Recent bison. This 
specimen was scored and the values obtained were 
plotted in the canonical variate space in Figure 12, 
to compare its position relative to that of Nearctic 
bison. The position of the east Siberian specimen, 
as one can see, is intermediate between athabascae 
(D = 2.85) and occidentalis (D = 3.23), suggesting 
a close affinity between these three forms. 

Crania of northern plains bison (CP, CGP, 
USRM, CRM) and athabascae were further com- 
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Figure 11. Projection of individual specimens and centroids of three taxa of modern bison onto the first two canonical variates, 
based on 11 cranial variables. Intercentroid distances (D) show the mutual relationship between taxa. Character vector lengths 
equal two standard deviations. 1 CHC; 2 GWA; 3 WHCO; 4 GPW; 5 AOW; 6 WMP; 7 RW2; 8 P2M3L; 9 M1M3L; 10 FP; 11 OC. 
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Figure 12. Three dimensional projection of the centroids of modern (1-6), early to mid-Holocene (7 and 8) and late Pleistocene 
bison (9) onto the first three canonical axes. The superimposed minimum spanning tree connects nearest neighbors. 1-athabascae; 
2-CP; 3-CGP; 4-USRM; 5-SGP; 6-EUS; 7-occidentalis; 8-antiquus. 9-priscus, The east Siberian bison (see text) was scored as 
an unknown to determine its affinities to the North American taxa. 


pared to those of subfossil Holocene bison from 
Alberta (Duffield, Fort Saskatchewan, Spirit 
River) using seven variables (Figure 14). The 
results, nearly identical to those in Figure 9, are 
interesting in that they show subfossil bison to be 
highly similar to modern athabascae with respect 
to the dimensions compared. 

Although none of the specimens used had been 
radiocarbon dated, other specimens from Duffield 
and Fort Saskatchewan that have been dated 
(McDonald 1981, Harington pers. comm.) suggest 
ages of between 4000 and 8000 BP for these sites. 

Harington (1984) described a cranial fragment 
from Banff, referred by him to B. b. occidentalis, 
which yielded a radiocarbon date of 3240+ 90 yrs 
BP. This demonstrates that occidentalis- 
athabascae-like bison lingered in western Alberta, 
an area occupied by B. b. bison in historic times, 
until relatively recently. 
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Figure 13. Dendrogram summarizing the distance relation- 
ships between modern (EUS, USRM, SGP, CP, CGP, 
athabascae), early to mid Holocene (occidentalis, antiquus), 
and late Pleistocene bison (priscus) in four-dimensional 
canonical variate space. The cophenetic correlation coefficient 
is 0.89. 
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Figure 14. Projection of individual specimens and centroids of geographic subsamples of B. bison and subfossil Holocene bison 
from Alberta onto the first two canonical variates, based on seven cranial characters, and showing minimum spanning tree with 


distances (D). 


To examine further the morphometric relation- 
ship of the modern forms of bison to earlier 
Holocene and late Pleistocene forms including 
those from Eurasia, a PCA on a correlation matrix 
of 8 cranial variables was done after standardiz- 
ing the data (Figure 15). The similar values and 
positive signs of the character loadings on the first 
PC axis show variation in size to be mainly along 
this axis. The second axis separates the forms on 
the basis of shape i.e. relative horn dimensions 
(SHTT, CLUC, TB) and cranial width measure- 
ments (GWA, WHCO and GPW). The Holocene 
and late Pleistocene representatives are situated to 
the right, while the smaller modern forms are 
predominantly on the left (athabascae is the 
notable exception). The forms with relatively 
wider skull dimensions and relatively shorter horns 
are found high along the second axis, those with 
longer horns and relatively narrower skulls lower 
down. The MST linking nearest neighbors, based 
on the distance matrix and projected onto the 
first two principal components, shows the cranio- 
metric relationships between the forms compared 
(Figure 15). 

The interesting points that emerge from this 
analysis are the divergent position of B. Db. 
athabascae among modern bison; its great similar- 
ity to Holocene forms and the general similarity 
of Holocene and late Pleistocene bison of Eurasia 
and North America. 
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The morphometric variation among modern, 
Holocene and late Pleistocene bison expressed in 
the above analyses provides a guide for the 
interpretation of the taxonomic relationships 
among them. The graphs in Figure 12 and 15 sug- 
gest that the variation within modern bison may 
be similar to that found among Holocene and late 
Pleistocene forms as well as that between modern 
and Holocene forms. For example the greatest dis- 
tance (D =7.68) between modern forms of bison 
in Figure 12, that between B. b. athabascae and 
EUS, exceeds that between occidentalis and anti- 
quus (D = 3.63) and that between occidentalis and 
priscus (D = 6.00); and is only slightly less than 
that between antiquus and priscus (D =7.80). In 
other words, variation within modern and late 
Pleistocene-Holocene bison is quite similar. If we 
now examine the differences between modern 
bison and late Pleistocene-Holocene bison 
(D 3.38-6.00), we find them to be of a similar 
magnitude to those found among modern bison 
(D 1.41-7.68). Similarly in Figure 15 the greatest 
distance (calculated following Goodman 1972) 
between modern forms is that between athabascae 
and bonasus (D=3.33), again exceeding that 
between occidentalis and antiquus (D= 1.56) as 
well as that between priscus and these two 
(D =2.37 and 3.19 respectively). The distances 
between geographic samples of modern bison 
range from 0.15 (CGP-CP) to 3.33 (athabascae- 
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Figure 15. Relationships among modern, Holocene and late Pleistocene bison from North America and Eurasia, on the basis of 
the first two principal components using eight cranial variables. The dots represent centroids of samples or individual specimens 
(Holocene or late Pleistocene specimen from eastern Siberia, Denmark, and the Netherlands) which are joined by a minimum 
spanning tree. Character vectors show the relationship of the eight original variables and the principal components. The first axis 
represents size variation, the second axis shape. 1 SHTT; 2 CLUC; 3 TB; 4 CHC; 5 GWA; 6 WHCO; 7 GPW; 8 ON. Data for 


Denmark and Europe centroid from Degerb@l and Iversen (1945). 


bonasus), those between modern and Holocene 
bison from 0.78 (EUS-subfossil European) to 3.71 
(bonasus-antiquus). Thus the degree of difference 
within the living subspecies is similar to that within 
late Pleistocene-Holocene bison as well as to that 
between modern and Holocene bison. 

If we take the amount of morphological differ- 
ence between the most divergent modern sub- 
species to indicate how much morphological 
difference may evolve without the acquisition of 
reproductive isolation (Mayr 1969), we are led to 
the conclusion that modern and late Pleistocene- 
Holocene bison are only subspecifically distinct. 

The proximity of the eastern Siberian specimen 
to occidentalis and athabascae in Figures 12 and 
15 is of particular interest, because of the sig- 
nificance of this specimen in interpreting the recent 
evolutionary history of North American bison. 
Univariate comparisons of 18 measurements of 
this specimen with those of the occidentalis and 
athabascae samples, using a ¢ test (Simpson et al. 
1960) showed none of the differences with 
occidentalis to be significant. Only two measure- 
ments differed significantly from those of atha- 
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bascae (WHCO, t=3.5335 P<.01 and WMP 
t=2.0192-P.< 305). 

The similarity between occidentalis and atha- 
bascae is further evident from the ratio diagram 
(Figure 16) comparing 19 cranial dimensions in 
B. b. bison, B. b. athabascae and B. b. occiden- 
talis. B. b. bison shows greater deviations from 
B. b. occidentalis than B. b. athabascae does. The 
differences are most marked in the horn measure- 
ments (SHTT, CLUC, TB, VD, TD, CHC) 
reflecting the general reduction in the size of the 
horns in B. b. bison. Of the other cranial dimen- 
sions the most obvious differences are exhibited 
in OC which is smaller in B. b. bison and in the 
tooth row measurements P2M3L, M1M3L which 
are relatively larger in B. b. bison. 

In B. b. athabascae the main difference from 
B. b. occidentalis is in the reduced horn length 
measurements (SHTT, CLUC, TB). Deviations 
from the occidentalis pattern are considerably less 
than in B. b. bison. The thickness of the horn core 
in B. b. athabascae (VD, TD, CHC) has changed 
relatively little compared to that of B. b. occi- 
dentalis, in contrast to the situation in B. b 


P2M3L - O S 


bison. The remaining cranial variables in B. b. 
athabascae are generally similar to those in 
occidentalis. 


4. Morphometric Analyses of Present Remnant 
Populations in Canada 


a) Craniometric comparisons of adult males. 
The taxonomic affinities of present Canadian 
plains bison herds (AP) relative to the historical 
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Figure 16. Ratio diagram comparing 19 skull dimensions in 
B. b. bison, B. b. athabascae and B. b. occidentalis, with 
the latter serving as a standard. 


bison populations were determined by scoring AP 
as an unknown in two CVA’s of historical popu- 
lations using the sample means of AP for 11 vari- 
ables (Figures 17 and 18). In both cases the AP 
mean values produced canonical variate scores 
close to the centroids of historical northern plains 
populations. Present Canadian populations of 
plains bison (AP) are therefore craniometrically 
representative of historical northern plains bison 
and close to the original population once found 
on the Canadian prairies. 

A CVA of adult AP and athabascae male refer- 
ence samples, the two ‘‘parental’’ populations of 
the present bison populations in Wood Buffalo 
National Park, was used to determine the taxo- 
nomic affinities of the NR and SE specimens by 
scoring them as ‘‘unknowns’’ and plotting their 
positions on the AP-athabascae axis (Figure 19, 
Table 2). In addition a three-group CVA was done 
including NR as an a priori group (Figure 20). In 
Figure 19 both NR and SE specimens occupy 
a somewhat intermediate position on the AP- 
athabascae axis. The NR centroid is closer to that 
of athabascae than to that of AP. Of the 16 adult 
NR males, 10 (62.5%) fall within two standard 
deviations of the athabascae centroid, only three 
(18.75%) fall within two standard deviations of 
the AP centroid and three (18.75%) fall in between 
the two standard deviation limits of the reference 
samples. In other words 81.25% of the NR speci- 
mens fall outside the two standard deviations of 
the AP centroid. 


Table 2 Canonical variate coefficients for 
11 cranial characters of adult male B. b. 
athabascae and AP bison (B. b. bison) 


Canonical variate 


Character coefficient 
1. ON -0.0658 
Zr VD 0.5633 
388 D 0.0730 
4. GWA 0.2993 
5. W -0.3765 
6. WHCO -0.0398 
T. GPW 0.1706 
8. P2M3L 0.0583 
9. MIM3L 0.4204 
10,.0€ 0.3536 
11. WMP -0.3254 


The scores for the five SE specimens present a 
similar picture. Two (40%) fall within two stand- 
ard deviations from the AP centroid, three (60%) 
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Figure 17. The morphometric relationship of present Canadian plains bison (AP) to historical populations of bison based on six 
horn and five cranial variables. Distances to nearest centroids are shown. 
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Figure 18. The morphometric relationship of present Canadian plains bison (AP) to historical populations of northern bison, 
based on nine cranial measurements and two horn measurements. Distances to nearest centroids are shown. 
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Figure 19. Histograms of canonical variate scores for four 
samples of adult male bison based on 11 cranial characters 
(see Table 2). AP and athabascae serve as reference samples. 
SE and NR were scored as unknowns on the AP-athabascae 
axis. The mean score for each sample with two standard 
deviations on either side is shown directly above the histo- 
grams. The distances of the NR centroid to the centroids of 
the reference samples is shown at the top of the graph. 
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fall within the two standard deviations range of 
the athabascae centroid. 

Three-group CVA produced quite similar results 
(Figure 20), with the NR sample being somewhat 
intermediate to AP and athabascae, but consider- 
ably closer to the latter. 

A number of cranial dimensions of adult 
athabascae, NR and AP males (Table 3) were also 
compared using a ratio diagram (Simpson et al. 
1960) with AP serving as a standard (Figure 21). 
The deviations of athabascae and NR from AP 
are on the whole remarkably similar with NR 
points paralleling those of athabascae, thus 
indicating that NR and athabascae skulls have 
nearly identical proportions and differ from AP 
skulls in similar ways. The position of most NR 
and athabascae points to the right of those of AP 
indicate larger size with athabascae generally being 
larger than NR. Differences in proportion lie in 
greater horn dimensions (SHTT, CLUC, TB, VD, 
TD and CHC showing the greatest deviations from 
AP) and relatively smaller WMP, CW, P2M3L, 
M1M3L and relatively larger OC. NR and atha- 
bascae differ significantly from AP in all measure- 
ments compared, except P2M3L and M1M3L (see 
Table 3). 

The single measurements that best discriminate 
between athabascae and AP skulls are FP, WHCO 
and VD with probabilities of misidentification of 
0.04 or less (Lubischew 1962) and corresponding 
percentages of joint nonoverlap of 96% or more. 
The best single discriminators between NR and AP 
are VD, CHC and FP with probabilities of mis- 
identification of approximately 0.05, 0.06 and 0.11 
and correspondingly lower percentages of joint 
nonoverlap than was the case in athabascae and 
AP. 

NR and athabascae differed significantly in 6 
out of the nineteen variables compared (OC, 
WHCO, GPW, FP, SHTT and CW) (Table 3). Of 
these OC, WHCO and GPW show the highest per- 
centage joint nonoverlap (72, 70 and 69% respec- 
tively). The proportion of the above variables in 
NR, although more like athabascae show a ten- 
dency toward the condition in AP. For example 
OC is relatively larger in NR than one would 
expect if they had been drawn from the same 
population as athabascae. The variation of the 
characters in the three samples is highly similar 
judging from the coefficients of variability in 
Table 3. None of the observed differences were 
significant. 
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Figure 20. Projection of individual specimens and centroids of AP, NR and athabascae males onto the first two canonical variates 
based on 11 cranial variates showing intercentroid distances (D). 


b) Craniometric comparisons of adult females. The ratio-diagram comparing 19 dimensions of 
Analyses of the adult female samples produced AP and NR adult female crania as well as those 
results similar to those obtained for the males. The of NMC 10405 (Table 5), shows deviations in NR 
two-group CVA of AP and NR specimens shows and athabascae from the standard (AP) similar to 
a similar degree of overlap to that found in the those in the male crania (Figure 24), but to a lesser 
males (Figure 22, Table 4). The generalized dis- extent. NR females differed from AP females sig- 
tance between NR and AP centroids is somewhat nificantly in 13 out of the 19 measurements com- 
less than that in the males. Of the 10 NR female pared in the ratio-diagram (Table 5). The best 
specimens 2 (20%) fall within two standard devi- univariate discriminators between NR and AP 
ations from the AP centroid, or 8 (80%) lie out- females were GPW, WHCO, RW, and GWA, but 


side the two standard deviations limits. The score 
for the athabascae female (NMC 10405) falls 


within one standard deviation to the left of the NR Table 4 Canonical variate coefficients for 
centroid and just outside the two standard devia- 11 cranial characters of adult female AP 
tion limit of the AP centroid. If the distance rela- (B. b. bison) and NR bison 

tionships of AP and athabascae cows were simi- 

lar to that observed in the males, the position of Canonical variate 
NMC 10405 would be toward the lower limit of Character coefficient 


ON -0.0904 
VD 0.4559 
TD -0.3054 
GWA 0.0482 


the athabascae distribution. The position of two 
SE and two WBNP females is similar to that of 
the SE males. 


i. 
2 
3. 
4. 

The three group CVA produced similar results 5. CW 0.0063 
(Figure 23) and underlines the existence of a mor- 6. WHCO 0.2168 
phometric relationship between NR and AP 7. GPW 0.0068 
females which is essentially the same as that found : a prion 
between the males of these groups. The position 10. OC 0.1089 
of the single athabascae female relative to the AP 11. WMP -0.2174 


and NR centroids is similar to that of the atha- 


bascae males in Figure 20. 
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the probabilities of misidentification were con- 
siderably higher than those found in the equiva- 
lent male samples (P~0.16, 0.18, 0.19 and 0.21 
respectively). Variation of the characters in AP 
and NR samples is quite similar (see coefficients 
of variation in Table 5). 

The single athabascae specimen closely parallels 
the NR sample, suggesting that NR female crania 
closely resemble those of the original B. b. atha- 
bascae. Comparison of the cranial measurements 
with the means for the NR sample showed that 
there were no significant differences (Table 5). On 
the other hand comparison with the AP sample 
demonstrated significant differences in five 
variates (GWA, WHCO, GPW, AOW, WMP, see 
Table 5). 
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Figure 21. Ratio diagram comparing 19 skull dimensions (see 
Table 3) in adult male athabascae, NR and AP with the latter 
serving as a standard. 
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NUMBER OF SPECIMENS 


c) Morphometric comparison of mandibles. 

CVA of 15 mandibular measurements of atha- 
bascae, AP and NR bulls yielded results similar 
to those for the cranial variables. The NR speci- 
mens are intermediate to AP and athabascae, but 
closer to the latter (Figure 25). Basic statistics for 
adult bulls and cows are listed in Tables 6 and 7. 


d) Sexual dimorphism of the cranium in AP and 
NR bison. 

Sexual dimorphism of the skulls of NR and AP 
specimens was analyzed by CVA of 11 cranial vari- 
ables. The degree of sexual dimorphism in the 
skulls of both groups was almost identical 
(Table 8). In both there is complete separation 
between male and female skulls, and the distance 
between the centroids of males and females is 
highly similar. The relative contribution to the dis- 


e WBNP 


OSE 
y, © ATHABASCAE 


-16 


CANONICAL VARIATE SCORE 


Figure 22. Histograms of canonical variate scores of three 
samples of adult female AP and NR bison based on 11 cranial 
variables. WBNP and SE specimens as well as the single 
athabascae specimen were scored as unknowns on the AP-NR 
axis. The mean scores for AP and NR samples, with two 
standard deviations on either side of the centroid is shown 
directly above the histograms. The distance (D) between AP 
and NR centroids is shown at the top of the graph. 
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Figure 23. Projection of individual specimens and centroids of AP, NR, athabascae, SE and WBNP females onto the first two 
canonical variates with intercentroid distances. The CVA was based on 11 cranial variables. SE, WBNP and the single athabascae 
specimen were combined as one a priori group, whose centroid is represented by the empty circle. 


Table 5 Means, standard deviations, coefficients of variation and observed ranges of 19 cranial 
variables in adult cows of NR and AP bison, with measurements of a single B. b. athabascae cow. 


Measurement athabascae NR AP t-test 


(N=1) N x SD CV_ Range N x SD CV _ Range NR-AP 
Basilar length (FP) 467 11 460.7 13.22 2.87 441-485 25 443 11.79 2.66 415-468 4.0006° 
Length, occipital crest to 217 15 208.6 9.68 4.72 188-226 33 209 8.66 4.14 184-225 n.s. 
nasal-frontal suture (ON) 
Spread of horns, tip-to-tip 435 14 459.1 37.61 8.19 401-515 25 438 40.89 9.33 359-517 2.63798 
(SHTT) 
Core length on upper curve, 147 15 146.2 27.84 19.03 103-182 34 123 23.22 18.90 63-162 3.04105 


tip-to-burr (CLUC) 

Length, tip of core to upper 134 15 132.2 20.37 15.41 98-159 34 115 19.69 17.19 63-145 2.78915 
base at burr (TB) 

Vertical diameter of core (VD) 51 16 50.1 5.98 11.93 43-64 35 44 4.98 11.30 30-52 3.8094¢ 


Transverse diameter of horn 54 16 51.9 6.27 12.07 41-62 a5 48 4.97 10.37 32-56 2.39264 
core (TD) 

Circumference of horn 164 16 163.1 18.06 11.07 139-199 34 148 22.76 15.35 100-248 2.33178 
core (CHC) 

Greatest width at auditory 220 14 218.3 9.87 4.52 205-240 24 205 6.98 3.40 189-220 4.8787¢ 


openings (GWA) 
Width of cranium between horn 227 16 224.6 9.20 4.09 210-242 34 209 8.10 3.88 190-222 6.1180¢ 
cores and orbits (WHCO) 


Greatest postorbital width 284 16 280.8 14.92 5.09 262-309 35° 257° YO.13 «393.94 237-275 6.8287°¢ 
(GPW) 
Anterior orbital width at 214 16 208.5 13.55 6.50 195-245 35 191 9.45 4.94 174-211 5.3527° 


notch (AOW) 
Width of skull at masseteric 177 16 170.4 8.50 4.99 155-186 36 «©6164 6.82 4.16 153-179 2.9285 
processes (WMP) 


Rostral width at maxillary- 114 14 101.5 6.91 6.80 89-112 36 99 5.64 5.67 83-110 n.s. 
premaxillary suture (RW1) 

Rostral width (RW2) 92 11 92.6 4.89 5.28 86-102 36 85 4.10 4.81 76-93 5.2054¢ 

Width of condyles (CW) 118 130 443585 G53) 25.73: “162-127 21 AN3*. 4 23-"3.76: 102-120 n.s. 

Width of occipital 128 15S 126.1 9.77 7.74 112-148 32 118 11.42 9.67 97-136 2.16954 
constriction (OC) 

P2M3 length (P2M3L) 141 16 138.9 11.34 8.14 123-154 36 147 4.49 3.06 137-158 3.7047¢ 

M1M3 length (M1M3L 89 16 92.2 5.49 5.95 86-109 36 92 3.36 3.64 84-99 fis. 

a P<.05 

oP = 01 

© P< .001 
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Figure 24. Ratiodiagram comparing 19 skull dimensions (see Table 5) in adult female athabascae, NR and AP bison, with the 
latter serving as a standard. 
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Figure 25. Projection of individual specimens and centroids of adult male AP, NR and athabascae onto the first two canonical 
variates with intercentroid distances (D). The analysis was based on 15 mandibular measurements. 


Table 7 Means, standard deviations, coefficients of variation and observed ranges of 15 mandibular 
variables in adult cows of B. b. athabascae, NR and AP bison. 


athabascae NR (N=S) AP (N= 37) 
(N:=.4) x SD CV Range X SD. CV Range t-test 
6 Length pm2-m3 (pm2m3L) 151 156.6 6.27 4.00 151-166 155.2 4.77 3.07 145-167 n.s. 
7 ml-m3 alveolar length 95 102.7 3.77 3.67 98-107 101.2 3.34 3.30 94-108 n.s. 
(mlm3L) 
9 Height of coronoid proc. 62 58.4 3.51 6.01 54-63 54.3 2.63 4.84 48-61 3.15106 
(HCP) 
10 Width of coronoid proc. 49 46.6 2.88 6.18 43-51 44.7 1.53 3.42 41-47 2.32004 
(WCP) 
11 Width of articulating proc. visi 78.6 4.72 6.00 72-83 74.0 3.01 4.06 69-81 2.9962 
(WAP) 
12 Posterior of m3 to mandi- 77 74.6 10.36 13.88 59-84 60.5 5.39 8.90 44-70 4.8728¢ 
bular foramen (Pm3F) 
13 Posterior of m3 to posterior 172 174.8 16.12 9.22 151-192 149.2 6.19 4.14 130-162 6.9092¢ 
edge of articulating proc. 
(Pm3AP) 
14 Width of ramus at posterior 28 29.6 1.67 5.64 28-32 28.6 1.72 6.01 26-34 n.s. 
of ml.(WR) 
16 Width of mandibular 48 49.2 4.38 8.90 43-55 47.0 2.08 4.42 42-50 n.s. 
condyle (WMC) 
17 m2 alveolar to base height | 60.4 2.30 3.80 58-64 59.8 3.26 5.45 54-68 n.s. 
(MD) 
18 Width of m3 (Wm3) 18 17.6 1.67 9.48 16-20 174 1.28 7.38 315-20 n.s. 
19 Width of m2 (Wm2) 18 18.0 2.00 11.11 16-21 18.0 1.30 7.22 15-20 n.s. 
20 Width of ml (Wm1l) 1g 17.8 0.84 4.71 17-19 17.0 0.95 ‘5.58: 15-19 n.s. 
21 Length of pm2-m1 (PL) 53 54.6 2.70 4.94 52-59 52.9 2.44 4.61 49-58 n.s. 
22 Substitute Total length (STL) 345 370.0 12.75 3.44 356-385 336.0 10.89 3.24 310-365 6.4342¢ 
4. P< .05 
bP 01 
SP 00! 
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2.51% 


crimination by individual variables differs some- 
what between the two samples on the basis of the 
standardized canonical variate coefficients. Over- 
all size and larger horn and cranial width dimen- 
sions in the males are the principal dimorphic 


Table 8 Standardized discriminant coefficients 
for 11 cranial variables in adult AP and NR 
bison and the distance (D) between the 
centroids for males and females, illustrate the 
nature and degree of sexual dimorphism of the 


pare chracteristics. 
The five best univariate discriminators (with 
AP NR approximate probabilities of misidentification) in 
ON 0.1378 -0.1642 AP bison are: GWA (P<.0003); WHCO 
vD 0.1796 0.4554 (P <.0007); VD (P=.003); GPW (P<.005); TD 
TD -0.0321 -0.5335 , 
GWA 0.4417 0.1703 (P~.006). In NR they are the same variables but 
Cw 0.0062 0.2100 in a somewhat different order WHCO (P< .003); 
WHCO 0.3786 0.3228 GWA (P< .008), VD (P< .01); GPW (P=.01) and 
GPW ~0.1529 0.3713 TD (P=.02). 
P2M3L 0.1395 -0.0784 
M1M3L -0.1892 -0.1432 e) Morphometric comparisons of the postcranial 
OC 0.0804 -0.3639 kelet 
WMP _ 0.0256 -0.0946 ee ay 
D oo-99 = 8.3 Doo-99 = 8.3 The plot of the first two principal components 
from the analysis of 25 postcranial variables shows 
500 
400 
300 
#*maie BISON 
% female BISON 
% co) } @maie NR 
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Figure 26. Relative position of specimens of B. b. bison, B. b. athabascae and NR in the projection on the first two principal 
components of a correlation matrix of 25 postcranial characters (see Table 8). Note the position of the single female B. b. bison. 
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a clear separation of plains bison on the one hand 
and athabascae and NR on the other into distinct 
clusters (Figure 26). The loadings on the first two 
components (Table 9) show those of the first com- 
ponent to be all of the same sign and nearly the 
same magnitude. We may, therefore, interpret the 
first component as predominantly a size compo- 
nent. The second principal component has load- 
ings of varying sizes and mixed signs, and can be 
said to summarize shape (Jolicoeur and Mosimann 
1960). The variation along the principal axes of 
the bivariate frequency ellipses represents size 
variation, with smaller individuals situated 
towards the lower left corner of the graph and the 
larger ones toward the upper right corner. Devia- 
tions at right angles of the principal axes repre- 
sent largely variation in shape. The different 
elevations of the plains bison and NR cluster 
relative to each other indicate differences in 
proportions of some of the postcranial elements. 
Examination of the loadings (Table 9) show that 
the elements of the front and hind limbs, particu- 
larly MCL, MTL and RADL, account for much 
of the difference in elevation. The lower position 
of the NR male sample indicates that NR bison 
have longer limbs than AP bison. The high posi- 


Table 9 Character loadings (correlation 
coefficients) of 25 postcranial variables on the 
first two principal components. 


I II 


1. SCAPL 0.984 0.118 
2. SCAPW 0.969 0.144 
3. HUML 0.998 -0.010 
4. HUMW 0.959 0.140 
5. RADL 0.979 -0.156 
6. RADW1 0.959 -0.124 
7. RADW2 0.988 -0.085 
8. MCL 0.955 -0.253 
9. MCWl 0.979 -0.098 
10. MCW2 0.984 -0.079 
11. FEML 0.988 -0.420 
12. FEMW1 0.984 -0.026 
13. FEMW2 0.988 -0.071 
14. TIBL 0.984 0.006 
15. TIBW1 0.988 -0.062 
16. TIBW2 0.964 -0.095 
17. MTL 0.955 -0.229 
18. MTW1 0.984 -0.055 
19. MTW2 0.984 -0.124 
20;-PELYV 0.979 0.063 
21. ILW 0.974 0.124 
22. ACETW 0.988 0.056 
23. ISCHW 0.935 0.210 
24. C7 0.833 0.449 
Zed | 0.935 0.282 


tive loadings of C7 and T1 indicate a vector for 
these variables in a direction approximately 
parallelling the major axes of the frequency 
ellipses. The position of the NR specimens rela- 
tive to the two specimens of athabascae suggests 
that the two are similar with respect to the post- 
cranial elements measured. 

Although only one plains female was available 
its position relative to the NR females suggests the 
differences observed in the males of both groups 
also hold true for the females. 

The univariate statistics for the postcranial 
variables are listed in Table 10. NR males differ 
significantly from athabascae males in only one 
measurement (width of the ischium) (Table 9). The 
plains bison bulls, on the other hand, differ 
significantly from those of athabascae in 16 and 
from those of NR in 22 of the 23 variables 
(Table 10). The lower number of significant 
differences with athabascae is probably a result 
of small sample size of athabascae. 

The best univariate discriminators as indicated 
by the probability of misidentification between NR 
and B. b. bison are MCL (P=0.05); RADL 
(P=0.09); RADWI1 (P=0.10); TIBL (P=0.10); 
HUML (P=0.12) and FEML (P=0.12) with esti- 
mated percentages of joint-non-overlap ranging 
from 88 to 95%. In athabascae and bison the best 
univariate discriminators are MCL (P=0.01); 
ISCHW (P=0.02); TIBL (P=0.05) and FEML 
(P~0.07) with percentages of joint-non-overlap 
ranging from 83 to 99%. 

Differences between NR and plains bison 
females (Table 11) were not tested because of small 
sample sizes. Examination of the means suggests 
however that the differences between the females 
are of the same magnitude as those observed 
between the males. 

The mean size of different postcranial elements 
in NR males is from 3 to 23% greater than those 
of NR females, with an average of 12.6%. In the 
plains bison the means for the same measurements 
in males exceed those of females by 5 to 22% with 
an average of 12.8%. Sexual dimorphism of size 
of postcranial elements in both groups is thus quite 
similar. 

Some of the postcranial variables were also 
compared in ratiodiagrams (Figures 27 and 28). 
The points for B. b. athabascae and NR bison 
show them to be highly similar to one another. 
Both differ from B. b. bison in size and in the 
proportions of certain postcranial elements, par- 
ticularly those of the pelvic girdle (e.g. ISCHW). 
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The lengths of the neural spines of the cervical 
(C7), and all thoracic (T) and lumbar vertebrae 
(L) of small samples of adult specimens in which 
the epiphyses had fully fused also showed signifi- ’ 
cant differences between the forms compared 
(Figures 29 and 30 and Table 12). The main differ- 
ences are found in the last cervical and the first Scary = 
twelve thoracic vertebrae. The neural spines of NR 
and athabascae specimens exceed those of bison 
and bonasus in length in diminishing amounts pane 
from front to back. In the male samples, NR and 
B. b. bison differ significantly in neural spine 
length of the first to 12th vertebra (Table 12). The FEMS 
difference is greatest in the second thoracic ver- 
tebra with an estimated joint-non-overlap of 92%. 
The neural spines of adult females are shorter than lite 
those of the males, but the differences between the 
female samples are similar to those observed in the 


NR 
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males, with NR and athabascae being very simi- hi 

lar to one another and greatly exceeding bison and icerw= 

bonasus in the length of the last cervical and the 

anterior thoracic spines (Figure 30, Table 13). acat 

Because of small sample sizes in the females the Figure 27. Ratiodiagram comparing 12 postcranial dimen- 
differences between NR, athabascae and bison sions of adult males NR, B. b. athabascae and B. b. bison, 
were not tested statistically. with the latter serving as a standard. 


Table 11 Means, standard deviations and observed ranges for 23 postcranial variables in adult 
cows of, NR and B. b. bison, and measurements of a single cow B. b. athabascae. 


B. b. athabascae(1) NR B. b. bison 
N X SD OR N x SD OR 
1. Scapula length 465 5 459.1 14.82 456-478 3 417.8 5.63 413-424 
2. Scapula width 232 4 224.8 5.87 219-233 3 205.3 9.29 195-213 
3. Humerus, length 342 5 350.8 13.98 336-369 3 327.8 16.38 313-345 
4. Humerus, width 44 5 43.3 Oe e 40-47 3 42.3 2.08 40-44 
5. Radius, length 314 5 332.1 13.69 324-356 Zz 308.0 22.62 292-324 
6. Radius, proximal width 91 5 91.1 3225 86-94 2 87.9 6.92 83-93 
7. Radius, distal width 80 5 82.1 4.60 77-86 2 80.7 5.30 77-84 
8. Metacarpal length 226 5 217.8 8.03 210-231 3 200.1 10.86 192-212 
9. Metacarpal proximal width 67 J 65.1 2.06 63-68 3 63.3 1.52 62-65 
10. Metacarpal distal width 65 5 Gi'.2 1.40 65-69 3 65.6 3.78 63-70 
11. Femur, length — 6 419.2 14.75 404-445 2 390.0 39.59 362-418 
12. Femur, proximal width — 6 135.3 3.79 131-142 Z 129.3 7.56 124-135 
13. Femur, distal width os 6 105.8 2.99 103-110 2 100.8 6.85 96-106 
14. Tibia, length 393 6 401.4 13.07 390-425 3 366.6 21.82 353-396 
15. Tibia, proximal width 109 6 111.3 3.42 107-117 3 110.6 6.65 103-115 
16. Tibia, distal width 70 6 70.7 1.84 68-73 3 70.0 1.0 69-71 
17. Metatarsal length 247 6 262.9 6.36 257-275 3 246.3 11.37 237-259 
18. Metatarsal proximal width 52 6 54.3 1.69 52-56 3 53.6 1352 52-55 
19. Metatarsal distal width 60 6 62.3 1:23 60-64 3 59.3 2.30 58-62 
20. Pelvis, length 484 5 477.3 11.67 458-487 1 416.0 _ a 
21. Width of ilium 461 4 446.2 21.84 422-474 2 429.4 3.67 427-432 
22. Width at acetabulum 253 6 247.4 7.60 235-254 pe 238.0 7.07 233-243 
23. Width at ischium 282 5 262.06 18.61 240-278 2 239.9 15.69 229-251 
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Figure 28. Ratiodiagram comparing 12 postcranial dimen- 
sions in adult female athabascae, NR and B. b. bison, with 


the latter serving as a standard. 
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f) Comparisons of external measurements and 
weights. 

Three external measurements (TL, T, SHHT) and 
body weight of athabascae, NR and AP animals 
of both sexes were compared by PCA after log 
transformation (Figure 31). The NR and AP dis- 
tributions show a relationship to one another simi- 
lar to that previously observed in the CVA of post- 
cranial variables. 

Here again the major axis represents the trend- 
line of increasing size and the difference in eleva- 
tion between NR and AP distributions indicates 
a difference in proportion, with NR specimens 
having greater TL, T and SHHT than AP animals 
of similar size. The measurements for the single 
athabascae specimen (the Radford bull) were 
derived from the literature (Radford 1911, Soper 
1941),! and may not have been taken precisely in 
the same way as in the others, but the proximity 
of this specimen to the NR specimens suggest their 
similarity. 


1 The length of tail cited by Soper (1941) appears to be erro- 
neous. An extract of Radford’s 1911 report in the Mammalogy 
Section files gives the length as 1’ 71/2’’. 
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Figure 29. Mean lengths of neural spines in adult male bison; B. b. bonasus (N = 9) (data from Roskosz 1962): B. b. bison (N= 4); 


NR (N=7) and B. b. athabascae (N=1). 


33 


500 


a= = 
._—- ooo 


400 uf , Ae BONASUS 


300 


200 


LENGTH OF NEURAL SPINES (mm) 


100 


C7, Fi ¥2 T3 T4 T5 T6 ™T7 #868 To T1O Tt Tia 13 -rae th L2 13 L4 L5 


VERTEBRAE 


Figure 30. Mean lengths of neural spines in adult female bison. B. b. bonasus (N = 13) (data from Roskosz 1962): B. b. bison 
(N=1); NR (N=2); B. b. athabascae (N= 1). 


Table 12 Greatest lengths of neural spines of the seventh cervical, the thoracic and the lumbar 
vertebrae of adult male bison. 


t-test 
B. b. athabascae NR B. b. bison NR v. B. b. bison 
N x SD OR N x SD OR N X SD OR 

C7 3S 420:7 6.35 418-429 7 354.7 24.26 307-383 4 391.0 34.18 352-434 n.s. 
Tl 4 547.0 13.66 533-565 7 520.4 22.36 483-544 4°7474.1- 20.17 457-501 3.41125 
T2 1” 542.0 7 532.1 21.96 491-552 4 474.1 18.55 453-493 4.43065 
v3 1 514 7 505.8 20.49 468-527 3 444.6 26.53 419-472 4.00305 
T4 1 488 7 47222 20.45 437-502 3 411.6 26.57 384-437 3.9666> 
TS 1 441 7 428.0 23.79 388-465 4 376.0 20.83 351-402 3.63146 
T6 i. 2992 7 381.8 22.14 345-417 3 2°335.5 2112 315-365 3.06515 
OF 1 350 7 340.1 18.61 310-367 3 303.0 19.14 286-330 2.8682? 
T8 1 311 7, 302s" * 13:20. 281-322 4 281.7 18.55 264-305 2.14204 
T9 D277 7 :272.4 10.2) 253-283 4 253.5 17.91 236-278 2.27024 
T10 242 WM ao 6.57 227-246 4 224.0 15.64 211-246 2.38484 
TI L. 210 7 204.2 4.82 196-210 4 188.7 10.43 179-202 3.4376> 
TZ i) a THEO 6.73 162-179 4 158.0 10.23 148-171 2.57104 
T13 L-~ 43 7 142.8 6.12 135-150 4° 134.7 8.80 126-145 mS. 
T14 borons2 t 12520 6.19 116-135 4-- 123.0 7.25 116-133 n.s. 
Ll be 130 7 Ass 7.17. 111-128 4 114.2 6.13 108-122 n.s. 
L2 1 ~126 7 108.0 6.16 101-120 4 = 102.7 6.13 94-107 Sy 
L3 tenis 7 101.0 6:55 95-114 4 95.5 5.06 88-99 nS: 
L4 1 — i) 94.2 7.65 87-110 4 89.2 5.67 81-94 n.Sz 
L5 l 95 7 86.2 7.31 81-100 4 _ 81.5 5.19 74-86 n.s. 
aP<.05 
% P<01 
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A similar comparison with two added measure- 95% confidence ellipse of AP males and are shown 


ments (RHT and HF) between AP and NR sam- to be absolutely larger in overall size and to have 
ples shows nearly identical results (Figure 32). NR relatively greater RHT, HF, TL and SHHT and 
adult males are situated largely (80%) outside the to a lesser extent T than AP males. 
0.3 
of % male AP 
& ~ @juv % female AP 
= ay - @ male NR 
it ra © female NR 
Pe % whee * Wy male ATHABASCAE 
as 


43 4.4 45 46 4,7 48 49 
0 
I 97.27 % 


Figure 31. Relative position of specimens of AP, NR and B. b. athabascae in the projection on the first two principal components 
of a variance-covariance matrix of four characters (logarithms of three external measurements and body weight). 


Table 13 Greatest length of neural spines of the seventh cervical, and the thoracic and the lumbar 
vertebrae of adult female bison. 


B. b. athabascae NR B. b. bison 


(N= 1) N X SD OR N X SD OR 
Ci 322 1 310 224 
Tl 422 2 415.5 7.77 410-421 2 375 1.41 374-376 
a2 421 1 420.0 1 379 
iE 400 2 400.0 2.82 398-402 ] 360 
T4 367 Z 374.5 0.71 374-375 1 344 
ES 333 2 343.5 Ps) 342-345 1] 321 
T6 303 Z 314.5 0.71 314-315 ] 296 
Tt 279 2 287.0 2.82.4, 285=289 ] 266 
T8 258 2 261.0 5.65 257-265 244 
T9 230 2 235.5 3.53 233-238 1 217 
T10 198 2 207.0 2.82 205-209 1 185 
Tll 163 2 b72.0 4.24 169-175 ] 155 
112 131 2 139-0 5.65 135-143 ] 127 
T13 110 2 118.5 4.94 115-122 107 
T14 95 Z 107.5 6.36 103-112 ] bi 
Ll 93 Zz O95 4.94 96-103 ] 90 
L2 85 1 90.0 i 82 
L3 80 2 74.5 0.70 74-75 1 75 
L4 74 2 70.0 1.4] 69-71 1 70 
E5 66 Z O75 Z.tZ 66-69 1 60 
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Figure 32. Relative position of centroids for different age classes of male and female AP and NR bison on the first two principal 
components of a variance-covariance matrix of six characters (Logarithms of five external measurements (1 TL; 2aT; 3 SHHT; 
4 RHT; 6 HF;) and body weight (WT)). Individual specimens as well as the centroid for adult males are shown for the NR sample. 
AP bison are represented by centroids of different age groups and the 95% confidence ellipses for males and females. The major 
axis of the ellipse represents variation in size, with large specimens in the lower right hand corner of the graph and small ones 
in the upper left. Centroids of the AP sub-samples are from left to right: two-year-old cows; three-year-old cows; adult cows; 
four-year-old cows: two-year-old bulls: four-year-old bulls; adult bulls; old bulls. NR bison are represented by from left to right: 
a three-year-old cow; a two-year-old bull; ten adult bulls. The character vectors show the relative magnitude and direction of 


the contribution of the six variables to the principal components. 


Comparison of the means for external measure- 
ments and weights show that NR and AP adults 
differ significantly in nearly all measurements 
(Table 14). NR bulls exceed AP bulls by 7 to 17% 
in the measurements compared and weigh on aver- 
age 17% heavier. Estimated joint-non-overlap for 
the measurements compared exceeds 80% except 
for tail length (Table 14). The cows differ from 
one another in a similar way with NR exceeding 
AP by 6 to 11% in the measurements compared. 
Estimated joint-non-overlap in the cows exceeds 
80% except for (Table 14). From the samples com- 
pared it appears that NR cows approximately 
equal AP bulls in some measurements. 

Measurements and weights of AP bison are 
similar to those of plains bison from Wichita 
Mountains Wildlife Refuge (Halloran 1960). Com- 
parison of T, HF, SHHT, and WT of adult AP 
bulls and 13 adult Wichita Mountains bulls (means 
T39.0; HF 59.96; SHHT 165.77; WT 589) showed 
significant differences in T (t = 3.6920 P<.01) and 
WT (t=5.0878, P< .001). The means for the same 
variables of 14 Wichita Mountains cows 
41 years and older, all differed significantly 
from those of AP. In all cases the means of the 
Wichita Mountains cows were smaller (T 34; 
HF 55; SHHT 139.9; WT 382). 


g) Comparison of outline shape and pelage 
characters. 
The mean differences in standardized outline 
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shape of adult plains, NR and European bison 
bulls and cows are shown in Figures 33 and 34. 
The most important differences in the profiles are 
in the lumbar area, the hump and the ventral 
contour of the head and neck. 

Both North American forms differ from B. b. 
bonasus in the more sloping outline of the pelvic 
and lumbar region, which is nearly level in bona- 
sus. The latter also shows a slight depression in 
the area of the neck in males (approx. 12 cm from 
the origin of the horizontal axis). NR bison differ 
from plains bison primarily in the shape of the 
hump, which reaches its highest point significantly 
farther forward in the former (Figure 33, 34). The 
highest point in the hump of the NR sample was 
10.35 cm from the origin, in the standardized out- 
line shapes compared to 8.85 cm in the plains male 
bison sample. The difference was highly signifi- 
cant (t= 6.7259 P< .001). A similar difference was 
found in the females with the mean for the highest 
point at 10.30 cm in NR and 8.82 cm in the plains 
females (t=2.1213 P<.05). The anterior outline 
of the hump in NR bison descends more precipi- 
tously than in the others, giving the hump a 
characteristic angular appearance. 

The differences in the humps of the three forms 
of bison are obviously related to the differences 
in the lengths of the neural spines discussed previ- 
ously. Other factors such as musculature, posture 
and differences in skeletal proportions (Mohr 
1939, Sokolov 1971) may contribute to the differ- 


Table 14 Means, standard deviations, coefficients of variation, and observed ranges for five 
external measurements (cm) and body weight (kg) in adult NR and AP bison. 


NR 66 AP 6d 
HB N 14 14 
X 305.5 253.0 
SD 13.81 19.92 
cy 4.52 7.87 
R 274-335 225-297 
T N 14 14 
Xx 49.6 44.3 
SD 3.20 4.19 
CV 6.45 9.45 
R 44-54 36-51 
Tail N 14 
length X 84.2 
incl., SD §.75 
tassel CV 6.82 
R 75-95 
HF N 14 14 
X 66.3 61.5 
SD 2.46 1.69 
CV 371 2.74 
R 62-71 58-64 
SHHT N 14 14 
x 182.2 165.5 
SD 9.30 7.58 
CV. 5.10 4.58 
R 168-201 152-179 
RHT N 14 14 
x 158.0 141.6 
SD 11.65 6.07 
CV 737 4.31 
R 144-181 137-152 
WT N 11 14 
+ 943.6 792.5 
SD 126.96 86.08 
CV 13.45 10.86 
R 759-1179 680-940 
aP<.0l 
b P<.001 


ences in outline shape. The conspicuous neck 
bulge of bonasus bulls, which gives many of them 
a double-humped appearance, (also seen in many 
bison depicted in European paleolithic cave art), 
appears to result from the underlying muscles 
(acromiotrapezius, splenius and rhomboideus) 
gaining greater prominence because of the shorter 
neural spines of the last cervical and the first two 
thoracic vertebrae in this taxon. In North Ameri- 
can bison the muscle bulge is obscured by the 
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t test NR 99 AP 99 t test 
8.10425 7 15 4.4179> 
254.5 231.3 
21.86 6.16 
8.58 2.66 
226-285 218-244 
3.7613> 7 15 n.s 
42.2 38.06 
3.81 4.78 
9.02 1235 
39-48 28-43 
8 
69.2 
10.33 
14.92 
55-91 
6.01755 6 15 4.5816 
61.5 ait aM | 
2.50 L335 
4.06 2.30 
59-66 56-60 
5.02315 iz 15 6.9443b 
162.2 146.6 
6.72 3.88 
4.14 2.64 
155-172 140-152 
4.67125 7 15 7.15185 
150.7 134.1 
aS 3.85 
4.74 2.86 
155-172 127-142 
3.54414 l 15 
508 497.6 
68.89 
13.84 
390-640 


greater development of the neural spines. Mohr 
(1939, p. 447) gives an alternative explanation for 
the double hump of B. b. bonasus attributing it 
to the head being carried higher and thereby push- 
ing the skin of the neck up to form a fold. The 
present analysis does not support this explanation. 
The differences in the dorsal outlines of the three 
bison are equally evident in both sexes, although 
the hump is more developed in the males. Aside 
from sexual dimorphism, age and to a lesser extent 


BONASUS 


Figure 33. Mean outline shape of adult bulls based on standardized outlines of B. b. bison (N= 17), NR (N= 14) and B. b. bonasus 
(N = 6). Only the dorsai and anterior ventral outlines (beard and neck mane) are shown as it is these areas that show the differences 


among the different forms. 


BONASUS 


Figure 34. Mean outline shape of adult cows based on standardized outlines of B. b. bison (N = 10); NR (N =9); and B. b. bonasus 


(N = 2). 


seasonal variation affect outline shape. Seasonal 
differences in outline shape of NR bulls consist 
chiefly in the reduction, and in some cases the 
virtual disappearance of the preputial tuft. 
Although seasonal changes in pelage are evident, 
loss of display hair similar to that described for 
B. b. bison by Lott (1979) was not observed. Age 
variation has however, a marked effect on outline 
shape. Old adults, particularly in the males, have 
more developed humps and heavier more robust 
bodies than younger adults (Plate 1). The reasons 
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for these differences lie in the greater length of the 
neural spines and greater muscular development 
of older bulls. The epiphyses of the neural spines 
are among the last to close. In B. b. bonasus 
epiphyseal fusion of the neural spines is reported 
to take place at age seven or shortly thereafter 
(Koch 1935). Although detailed studies of the 
chronology of epiphyseal fusion in North Ameri- 
can bison are still lacking, it appears highly prob- 
able that it is similar to that in B. b. bonasus 
(Duffield 1973). 
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Ventrally, the outline shows the greater develop- 
ment of the beard in the North American bison, 
(particularly in B. b. bison) compared to B. b. 
bonasus. The neck mane in NR bison is generally 
less developed than is the case in B. b. bison and 
B. b. bonasus. The general effect of this is that 
the neck in adult NR bison tends to be more dis- 
tinct. The ventral outline is thus determined by 
differences in hair length, which brings us to the 
analysis of pelage characters. 

After preliminary exploratory evaluation of 
pelage characters in AP and NR bison, the number 
of characters selected for further analysis was 
reduced to four. The criterion for selection was 
the presence of two or more recognisable and 
definable character states that lent themselves to 
reasonably objective evaluation by an observer in 
the field or from good photographic images. 

The four pelage characters selected and their 
states and their estimated frequency of occur- 
rence in AP and NR populations are shown in 
Figure 35. Three of the four pelage characters 
(2, 3 and 4) show clear polarity of the character 
states between NR and B. b. bison, while the cape 
shows a greater degree of overlap. 

My observations on pelage characteristics and 
dorsal outline of NR and AP bison generally agree 
with those of Geist and Karsten (1977) and Rey- 
nolds et al. (1982), but the distinctions listed by 
them are not always as clear cut as they lead us 
to believe and there are some discrepancies. My 
evaluation of some pelage characters listed by 
them was not conclusive because of inadequate 
data. 

The differences in colour of the lower neck and 
the cape were not found to be as distinct as those 
authors suggest. Colour appeared to be quite vari- 
able in both NR and AP specimens and there was 
considerable overlap. However, light (blond) capes 
showing a clear cut demarcation appear to be 
characteristic of plains bison. 

According to Geist and Karsten (1977) the 
preputial tuft of NR bulls appears to be shorter 
and thinner than those of AP bulls. I was not able 
to confirm this because I lacked measurements for 
AP bulls. Similarly tail length including the tassel 
could not be compared to evaluate Geist and 
Karsten’s (1977) impression that the tail of NR 
bison was longer and hairier as no data were avail- 
able for AP specimens. My observations indicate 
that these differences are not very great. Tail 
length is much more variable than other external 
measurements (Table 14). In both NR and AP 
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bison the end of the tassel may reach the hock. 
However, in some prairie bison herds, presuma- 
bly those with a greater admixture from southern 
stock, the tail may be noticeably reduced in length. 
NR and plains bison according to Geist and 
Karsten (1977) differ in that the contour of the 
back from the top of the hump to the root of the 
tail forms a descending line in the latter but not 
in the former. Outline shape comparison in the 
present study (Figure 33 and 34) showed NR and 
AP specimens to have similar descending dorsal 
outlines, but significant differences in the outline 
of the hump. Similarly, the outline shape compar- 
ison did not show any difference in the length of 
the neck of NR bison as suggested by Reynolds 
et al. (1982). 

Geist and Karsten (1977) thought NR bison to 
be intermediate in pelage characteristics between 
wisent and prairie bison and postulated that juve- 
nile pelage of NR bison would be intermediate in 
its characteristics between that of wisent and plains 
bison. One NR calf photographed on June 30, had 
a reddish juvenile pelage similar to that of the 
prairie bison and unlike the greyish brown juve- 
nile pelage of the European bison. Some other 
calves observed in the latter half of August were 
largely dark but showed remnants of the reddish 
juvenile pelage. It appears, therefore, that NR 
bison resemble the prairie bison with respect to this 
character. Available data did not allow a detailed 
comparison of the juvenile pelages of NR and 
plains bison. 

The same authors state that NR cows are less 
sexually dimorphic than plains cows, being closer 
to bulls in body size and proportions, as well as 
in horn form and pelage characteristics. They 
found NR cows on the whole more difficult to 
identify as cows. The analyses in the present study, 
however, have demonstrated that there is no 
appreciable difference in the degree of sexual 
dimorphism in NR and AP bison with respect to 
body size and proportions or horn form. As to 
pelage characteristics there is possibly a more pro- 
nounced sexual dimorphism in the development 
of display hair in AP bison. The greater difficulty 
in distinguishing the cows from bulls in NR bison 
is related chiefly to the female’s greater overall size 
and high hump. As we have seen the body size- 
shape relationship in male and female bison is 
essentially the same, with centroids along a com- 
mon axis (Figure 32). It follows from this that 
large cows tend to have a more bull-like appear- 
ance than small cows. 


The differences in outline shape and pelage 
characters observed in NR bison in all probability 
reflect the character states typical of the original 
athabascae. Unfortunately, there is little informa- 
tion on the external appearance of the original 
athabascae, thus precluding a detailed comparison 
of the external pehnotypes of NR and athabascae. 
However, what little information does exist, con- 
firms the above conclusions. The two photographs 
of a wood bison bull taken by Kitto (Graham 
1923, Kitto 1924) show the characteristic hump 
outline seen in NR bulls. Seton (1927) commented 
on the absence of a mane, beard and ‘‘chaps’’ in 
the Radford bull (collected in December). The type 
specimen (NMC 299) (collected in March) likewise 
shows little development of the mane and beard 
and has no ‘“‘chaps’’, and a sketch by Rowan 
(University of Alberta Archives), who collected 
two wood bison bulls in Wood Buffalo National 
Park in October 1925, shows similar character 
states. As two of the above specimens were col- 
lected in winter, it could be argued that wear 
caused the observed states. However, NR bison 
observed in Elk Island National Park, did not 
show great seasonal differences in outline shape. 
Seasonal differences consisted primarily of a 
lengthening of body hair and a conspicuous reduc- 
tion in the prominence of the preputial tuft. It is 
likely that the effects of abrasion on the pelage 
would be most noticeable in late winter. 

As to colour, the original athabascae is gener- 
ally described as darker than the prairie bison 
(Seton 1886, McFarlane 1905, Graham 1923, 
Seibert 1925), but Mignault (Seton 1886) states 
that in summer the wood bison assumed a hue 
similar to that of the prairie bison. 

Although, we may thus regard the present NR 
bison as close to original athabascae, we could 
expect certain plains bison traits to be evident on 
the basis of the previous osteological evidence. The 
frequency of occurrence of typical plains character 
states in the NR sample analysed appear to bear 


this out (Figure 35). The presence of chaps and 
plains bison type frontal display hair are of 
particular interest. Although hybridization exper- 
iments between North American and European 
bison and various other bovini have been con- 
ducted since the 1920’s (Gray 1972), detailed infor- 
mation on the inheritance of certain traits is 
lacking. We do know, however, that in prairie 
bison-wisent F, hybrids, the prairie bison traits of 
chaps and frontal display hair dominate (Mohr 
1952). We may assume that holds true also for 
plains-wood bison crosses. If these traits are deter- 
mined by a single pair of alleles, as seems probable 
(Deakin et al. 1935), the frequency of the plains 
allele in the NR population, based on the esti- 
mated frequency of occurrence of the phenotypes 
in the population and the Hardy-Weinberg for- 
mula, would be about 0.05. If the above assump- 
tions are correct, the allelic frequency would indi- 
cate the extent of possible introgression of the NR 
population by AP genes. 

A weighted character index using four pelage 
characters and one outline shape character (the 
hump) was devised (Appendix 3) as a means of 
evaluating the external phenotypes in the field or 
from adequate photographic records. A sample of 
33 AP and 32 NR adult bulls were scored. Of the 
AP sample 18 (54%) had scores from 90 - 100. 
13 (39%) from 80 - 90 and 2 (6%) from 70 - 80. 
For NR bulls the scores were as follows: 9 (28%) 
had scores from 0 - 10; 20 (62.5%) from 10 - 20; 
one (3%) from 30 - 40 and two (6%) from 50 - 
60. NR bulls with high scores are characterized by 
plains characteristics (Plate 1). The greater range 
of scores (0 - 60) in the NR sample indicates a high 
degree of variability possibly a result of introgres- 
sion of B. b. bison genes. 

Examples of variation in external characters 
observed in male and female NR bison are illus- 
trated and compared to typical AP bison in 
Plates 1 to 4. 


Plate 1. Photographs illustrating the range of phenotypic variation in NR bison bulls. 


1. Fully adult NR bull showing a combination of external characters 
considered to be typical of B. b. athabascae. (character index 
score 0). Note the robust heavy body characteristic of fully mature 
males, the dark colour, procumbent and thin frontal display hair, 
thin beard, little developed ventral neck mane, lack of clear demarca- 
tion of the cape, lack of chaps, robust horns with slightly spirally 
twisted tips and dorsal outline with highest point far forward. Photo 
taken in Elk Island National Park, 28 July 1982, by Donna 
Naughton. 


2. Full grown adult NR bull (character index score 10) showing a 
combination of characters similar to that of the previous bull, but 
having a somewhat lighter cape. Photo taken in Elk Island National 
Park, 3 October 1984, by Donna Naughton. 


3. A younger NR bull (character index score 20) characterized by 
a less robust body and lower neural vertebral spines. Pelage charac- 
ters and outline shape are those of B. b. athabascae. Photo taken 
in Elk Island National Park, August 1977, by David B. Campbell. 


4. An even younger NR bull already showing the characteristic 
outline shape and procumbent frontal display hair. Photo taken in 
Elk Island National Park, 5 October 1984, by Donna Naughton. 
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Plate 1. Photographs illustrating the range of phenotypic variation in NR bison bulls. 


6. Adult NR bull showing some plains bison characteristics (chaps, 
relatively well developed frontal display hair, hump more like plains 
bison). Photo taken in the Metro Toronto Zoo, 3 September 1984, 
by C.G. van Zyll de Jong. 


8. Adult male plains bison (AP) (character index score 100) showing 
the typical northern plains bison phenotype. Note dense standing 
frontal display hair, well developed ventral neck mane, chaps, light 
cape with clear posterior demarcation contrasting with dark head, 
dorsal outline (‘‘carp back’’) and relatively thin simple horn. Photo 
taken in Waterton Lakes National Park, in late summer by 
J.C. Holroyd. 


5. Adult NR bull (character index score 40) showing a combination 
of plains (frontal display hair, chaps) and wood bison (hump) charac- 
teristics. Photo taken in the Moose Jaw Wild Animal Park, 17 June 
1983, by H.W. Reynolds. 
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7. A younger NR bull from the same herd showing similar traits 
as well as a blond cape. Photo taken in the Metro Toronto Zoo, 
3 September 1984, by C.G. van Zyll de Jong. 
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Plate 2. Phenotypic variation in NR bison bulls 


.. Profile of adult NR bull (same individual as in Pl. 1, #1) showing 2. A close up of the head of the same individual. 
ypical procumbent frontal display hair, thin pointed beard, poorly 

leveloped neck mane and uniform dark cape. This individual in all 

ikelihood embodies the pure primitive phenotype of the original 

3. b. athabascae. Photo taken in Elk Island National Park, 28 July 

982, by Donna Naughton. 


. Adult NR bull, a somewhat smaller and probably younger animal 4. Close up of the same individual. 
han the bull in 1 and 2 but showing the same phenotype. Photo 

aken in Elk Island National Park, 28 July 1982, by Donna 

Jaughton. 


5. Adult NR bull with distinct plains bison phenotype 
(character index score 60). Note the dense erect frontal 
display hair, blond cape clearly set off from the rest of 
the body and simply curved thin horn. Photo taken in 
Wood Buffalo National Park, 10 September 1967, by 
E. Kuyt. 
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Plate 2. Phenotypic variation in NR bison bulls 
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6. Head of an adult NR bull having phenotype similar to the animal 7. Front view of the same NR bull as in fig. 3 and 4, representing 
in the previous picture. Photo taken in Elk Island National Park, the athabascae phenotype (character index score 0). 
5 October 1984, by Donna Naughton. 
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8. Front view of adult NR bull (character score 20) with somewhat 9. Front view of NR bull (character index score 40) with similar cape 
more erect frontal display hair and lighter cape. Photo taken in Elk and frontal display hair confirmation as the bull in figure 8, 
Island National Parks, August 1977, by David B. Campbell. but with heavy chaps. Photo taken in the Metro Toronto Zoo, 


3 September 1984, by C.G. van Zyll de Jong. 
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Plate 3. Phenotypic variation in NR bison cows compared to AP cows. 


1. Adult NR cow representing athabascae phenotype (character index 
score 0). Note lack of chaps, short sparse hair on head and dorsal 
outline. Photo taken in Elk Island National Park, August 1982, by 
Donna Naughton. 
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2. Adult NR cow showing similar phenotype. Photo taken in Elk 
Island National Park, 5 October 1984, by Donna Naughton. 
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3. Adult NR cow with relatively well developed beard and ventral 
neck mane. Photo taken in Elk Island National Park, 5 October 
1984, by Donna Naughton. ) 


4. Adult NR cow with light cape and chaps. Photo taken in the Metro 
Toronto Zoo, 3 September 1984, by C.G. van Zyll de Jong. 
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Plate 3. Phenotypic variation in NR bison cows compared to AP cows. 


5. NR cow with the characteristic forward athabascae hump. 
Photo taken in Elk Island National Park, October 5, 1984, by 
H.W. Reynolds. 


6. NR heifer already showing characteristic forward hump. Photo 
taken in Elk Island National Park, 5 October 1984, by Jo-Anne 
Reynolds. 


7. Adult plains bison cow (AP) showing characteristic ‘‘hair-roll’’, 
well developed beard and ventral neck mane, chaps and ‘‘carp’’ back. 
Photo taken in Elk Island National Park, 5 October 1984, by Donna 
Naughton. 


& 
P 
8. Adult plains bison cow (AP) showing similar characters. Photo 
taken in Elk Island National Park, 5 October 1984, by Donna 
Naughton. 
¥"* 
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1. Head of adult N 

fig. 1). Showing short hair on the head and poorly devel- 
oped beard and neck mane, characteristic of the athabascae 
phenotype. 
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3. Adult NR cow showing moderate amount of frontal hair 
and clean legs. Photo taken in Elk Island National Park, 


August 1977, by David B. Campbell. 


Plate 4. Phenotypic variation in NR bison cows compared to plains bison c 


ows (AP). 
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2. Old adult NR cow, note procumbent frontal hair, short 
beard, clean legs and long right horn (left horn is broken). 
Photo taken in Elk Island National Park, 28 July 1982, by 
Donna Naughton. 
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4. Face of adult NR cow showing short frontal hair and 
beard. Photo taken in the Metro Toronto Zoo, 3 September 
1984, by C.G. van Zyll de Jong. 


Plate 4. Phenotypic variation in NR bison cows compared to plains bison cows (AP). 
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5. Adult NR cow with relatively large amount of pro- 
cumbent frontal hair, small beard, and forward hump. Photo 
taken in Elk Island National Park, 5 October 1984, by | % 
H.W. Reynolds. , =” 


6. Face of adult NR cow showing similar characteristics. 
Photo taken in Elk Island National Park, August 1977, by 
David B. Campbell. 


7. Adult plains bison cow (AP) showing a well developed 8. Adult plains bison cow (AP). Note the well developed 
“hair roll’’ characteristic of most plains bison cows. Photo frontal hair, beard and chaps. Photo taken in Elk Island 
taken in Elk Island National Park, 5 October 1984, by Donna National Park, August 1977, by David B. Campbell. 
Naughton. 
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Conclusions and Discussion 


1. Historical Geographic Variation and 
Subspeciation of North American Bison. 


The preceding craniometric analyses of samples 
representing historical North American bison 
populations show that geographic variation was 
largely clinal in a predominantly south — north 
direction. However, a phenotypic discontinuity is 
revealed between boreal populations (B. b. 
athabascae) and plains populations (B. b. bison), 
coinciding approximately with the ecotone 
between grasslands and boreal forest in the north- 
western part of the historical distributional range 
(Figure 10). This discontinuity, involving size as 
well as shape of the skull, fully justifies the recog- 
nition of the subspecies B. b. bison and B. b. 
athabascae. Evidence from analyses of postcranial 
morphology and external characters further sup- 
ports the phenotypic discontinuity. 


McDonald (1981) interpreted the differences 
between the two subspecies to be largely those of 
size, aside from some minor difference in horn 
core shape, and the transition between the subspe- 
cies merely as an increase in the slope of the north 
— south character size cline. However, of the 
change in six cranial variables with latitude shown 
graphically by him only two, spread of horn cores 
and horn core length along the upper curve, show 
a definite clinal change. The other four are prob- 
ably not clinal. The change in limb characters 
(McDonald 1981), including humerus, radius and 
metacarpal length, are even less convincingly 
clinal. To the contrary, all show relatively little 
change between samples of plains bison and a con- 
siderable change between the northern plains 
population and B. b. athabascae suggesting a 
discontinuity rather than a cline. 

Specimens from the Canadian Rocky Moun- 
tains, as well as those from areas west of the great 
plains in the United States are morphometrically 
close to those from the prairies. It is highly prob- 
able that prairie bison penetrated and crossed the 
mountain barrier through the predominantly east- 
west oriented river valleys and mountain passes. 
The montane and intermontane bison populations 
are thus considered to be B. b. bison. This con- 
clusion agrees with that of McDonald (1981) and 
is contrary to that of Skinner and Kaisen (1947). 
The last mentioned authors included those popu- 
lations in B. b. athabascae on the basis of larger 
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measurements, and they postulated a dispersal 
route from the north through valleys between the 
Rocky Mountain ranges. 

Cranial material from the southern Peace River 
region in Alberta and British Columbia analysed 
in this study, yielded somewhat ambiguous evi- 
dence as to the affinity of the historical bison 
population in that area. The analysis of six horn 
core and five cranial measurements of a small 
sample indicates near identity with the historical 
population on the Canadian prairies (Figure 6). 
The analysis of a different set of characters, com- 
prising five cranial and two horn dimensions, 
resulted in a more intermediate position of the PR 
centroid, between plains bison and wood bison, 
although still considerably closer to the former 
(Figure 9). Finally, analysis of a much larger 
sample, but using only three variables, shifted the 
PR centroid closer to athabascae and away from 
bison (Figure 10). How are we to interpret these 
seemingly contradictory results? In my opinion the 
most likely explanation is that the population in 
the upper Peace River area was somewhat inter- 
mediate between B. b. bison and B. b. athabascae, 
but resembling the latter taxon more closely. The 
first two analyses, based on only five specimens, 
may not have adequately reflected variation in the 
upper Peace River populations. The larger sam- 
ple used in the analysis of three variables, on the 
other hand, probably represents a more realistic 
picture. The area between the upper Peace River 
parklands and the prairies proper, consisting of 
heavily wooded country, constituted the transition 
zone between the two subspecies. 

This conclusion is further supported by the 
information given to Seton (1886) by Elzear 
Mignault who lived in the Peace River area from 
1863-75. According to Mignault the bison in the 
Peace River area were distinct from the prairie 
bison. He also implied that prairie bison occasion- 
ally penetrated the area, the last being killed in 
1866. The area where the ‘‘prairie buffalo’’ was 
killed is not specified, although it is likely the 
observation was made in the valley of the upper 
Peace River. 

Historical references to so called wood buffalo 
in Manitoba (Roe 1970) are not backed up by evi- 
dence from specimens. Here, as well as in other 
instances (see Roe 1970), we seem to be dealing 
with the uncritical application of the term ‘‘wood’”’ 


to plains bison observed in wooded country. 
The historical distribution of B. b. athabascae 
was thus, on the basis of the objective evidence, 
restricted to the northwestern portion of the spe- 
cies’ range (Figure 36), an area closely coinciding 
with that outlined by Soper (1941). The exact 
boundaries of the historical distributional range 
of B. b. athabascae cannot be determined at 
present on the basis of the available information. 
The distribution of specimen localities, localities 
given in historical reports and the distribution 
of toponyms relating to bison suggest the pos- 
sible existence of a zone between the two taxa 
(Figure 36) that may have been only intermittently 
or sparsely populated. | 
Subfossil material referred to B. b. athabascae 
indicates a prehistoric range extending further 
north and west into the Yukon and Alaska 
(Figure 36). Three radiocarbon dates of 1350, 
1465 and 1810 yrs BP for B. b. athabascae have 
been reported for the Yukon (Harington 1980a). 
Two localities from which B. b. athabascae 
specimens have been reported, from near Fort 
McPherson on the Peel River (Skinner and Kaisen 
1947, MacDonald 1981) and from western British 
Columbia (Smith 1977) are here rejected on 
grounds of morphometric and other evidence. 
The Peel River locality is based on a partial 
cranium of an adult bull (USNM P 16861) col- 
lected by E.A. Preble in 1904. The discriminant 
score for this specimen (90.2) based on nine vari- 
ables (Appendix 2) identifies it as a B. b. bison. 
The probability that the specimen was drawn from 
the same population as the athabascae reference 
sample is exceedingly small (P< .0001). A note by 
Preble attached to the specimen and not com- 
mented on by the authors previously named fur- 
ther undermines confidence in the origin of this 
specimen. The note reads as follows: 


“*Bison skull obtained at Fort Simpson, Mackenzie, 1904 
by Edward A. Preble. 

“‘This specimen was in one of the warehouses of the 
Hudson’s Bay Company at Fort Simpson, having been 
saved probably because a sort of museum was for 
many years maintained there. The skull was stated by 
A.F. Camsell who had lived at that post for a number 
of years, to have been brought from Peel River (Fort 
McPherson). However in view of the fact that it bore 
no label and might have been confused in Mr. Camsell’s 
mind with some other specimen it is quite possible that 
it really came from somewhere about Great Slave Lake, 
where the bison was common to within a comparatively 
recent date, and near which it still occurs. [signed] 
EA-P.’’. 
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It seems highly probable that the skull was 
brought in from a locality much farther south than 
Preble suggested. 

The skull from coastal British Columbia 
reported on by Smith (1977) was also identified 
by me as a B. b. bison on the basis of its published 
measurements. The unusual locality of the find, 
near Kitimat in the coast forest, moreover, casts 
doubt on the specimen’s origin. 

The proposed recognition of two subspecies in 
the plains bison (Krumbiegel 1980) is not sup- 
ported by the results of the craniometric analyses. 
Reliable information on geographic variation in 
external characters of historical plains bison popu- 
lations is sparse. Aside from some old photo- 
graphic records, there appear to be no data that 
can be used in forming an idea of geographic and 
individual variation of external characters in 
historical prairie bison. A few old illustrations can- 
not be considered as an adequate foundation for 
the designation of a subspecies, as Krumbiegel 
(1980) has done. Moreover, it is reasonable to 
assume, given the continuous range across the 
Great Plains, that variation in pelage characters, 
like that of skeletal characters was clinal in nature. 
No useful purpose is served by taxonomically sub- 
dividing contiguous populations showing clinal 
variation. 


2. Taxonomic Status of Extant North American 
and European Bison 


The multivariate craniometric comparison of the 
two North American subspecies and the European 
bison lend support to the view that all living 
bison are conspecific (Bohlken 1967). The closer 
craniometric proximity of B. b. bison and B. b. 
bonasus to one another, than that of either to 
B. b. athabascae is difficult to harmonize with the 
view that the European bison is specifically dis- 
tinct. The results of the craniometric analysis 
leaves us two options — either there are three 
species present or the scatter of the three groups 
represent intraspecific variation. Independent evi- 
dence from ecological and behavioural studies 
(e.g. Fuller 1960, Heptner et al. 1966, Krischke 
1984, Lott 1974, 1979, McHugh 1958, Sokolov 
1979) of American and European bison as well as 
the complete interfertility of the two favour the 
latter option. If one considers all available evi- 
dence, and one applies the criteria proposed by 
Mayr (1969 p. 197) to the question ‘‘subspecies or 
allopatric species?’’ one is led to conclude that 


Figure 36. Historical distribution of the bison in northwestern North America based on specimen localities (S), historical reports 
(R), and toponyms (T), and showing the probable subspecific boundary between B. b. bison and B. b. athabascae. Localities 


where prehistoric bison referred to B. b. athabascae (F) have been found indicate the possible extent of the prehistoric distribution 
of this taxon. 
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American and European bison are subspecies of 
one species. Bohlken (1967), who compared the 
European bison and the American plains bison by 
means of univariate and bivariate techniques 
arrived at the same conclusion. I fully concur with 
his argumentation and taxonomic conclusions. 
Reluctance to accept conspecific status of the 
European and American bison is undoubtedly 
influenced by the marked phenotypic differences 
between them and their extreme disjunct histori- 
cal distributions. These are, however, of relatively 
recent origin. The phenotypic differences dis- 
played by modern forms of bison are, moreover, 
no greater than those observed in some other 
Holarctic species, e.g. Cervus elaphus, with simi- 
lar late Pleistocene-Holocene distributions, and 
less than those found between the most divergent 
subspecies of Syncerus caffer (Grubb 1972). These 
examples illustrate the possible extent of variation 
within a species and show that the degree of differ- 
ence among Recent bison is within the observed 
limits of intraspecific variation of related species. 


3. Morphometric Relationships of Late 
Pleistocene-Holocene Bison and Subspecific 
Differentiation of North American Bison 

The morphometric analysis of bison crania from 
the late Pleistocene and Holocene sheds some light 
on the possible origin and relationships of living 
forms. The comparisons show the well known 
progressive reduction in size (Holocene dwarfing 
Wilson 1980), as well as variation in the shape of 
the cranium and horn cores. An interesting obser- 
vation that emerged from these analyses is that the 
variation within modern bison may be similar to 
that found within and between Holocene and late 
Pleistocene forms, as well as that between modern 
and earlier forms of bison. Viewed together in 
multivariate space they form one morphometric 
continuum, which is best regarded as representing 
one variable chronospecies. Thus late Pleistocene 
and Holocene bison from different parts of 
Eurasia were not unlike North American occiden- 
talis, and the derivation of modern bison from a 
common stock as suggested by Skinner and Kaisen 
(1947), Fuller and Bayrock (1965) and Guthrie 
(1970, 1980) appears highly likely. 

The presence of athabascae-like bison in east- 
ern Siberia at the close of the last glaciation 
(Flerov pers. comm.) is particularly relevant to the 
question of the subsequent evolution of North 
American bison. This small-horned eastern 
Siberian bison is represented in the present analysis 
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by measurements from a well-preserved skull from 
the Bol’shaya Chukochya River in the Kolyma 
River lowlands (Sher 1971). The similarity of this 
bison to North American occidentalis and to 
athabascae is strong evidence of the former exis- 
tence of an ancestral Beringian population. After 
the reestablishment of the Bering Strait at the 
beginning of the Holocene, the Siberian popula- 
tion became extinct, probably during the Hypsi- 
thermal (Flerov 1979), but the American popula- 
tion survived and expanded southward. The 
disappearance of the ice barrier eventually led to 
contact and hybridization with the southern form 
(B. b. antiquus), as intergrades between the two 
indicate (Wilson 1975, Guthrie 1980, Harington 
1984). The modern prairie bison (B. b. bison) 
probably evolved from this mixed stock during the 
Hypsithermal (about 7000-5000 BP) when exten- 
sive grasslands appeared (Harington 1984), and 
gradually increased its range as the grasslands 
expanded. At the same time woodland and park- 
land adapted forms (occidentalis-athabascae) were 
increasingly restricted to the northwestern parts 
of the species’ range. Subsequent climatic, and 
possibly other factors as well, eventually resulted 
in contraction of athabascae range to its histori- 
cal proportions (Figure 36). B. b. athabascae may 
be considered the more primitive of the two North 
American subspecies on the basis of its greater 
morphometric similarity to B. b. occidentalis and 
the retention of primitive external characters. 
Seton (1886) was the first to express the view that 
athabascae was the more primitive form, a view 
subsequently expressed or implied by Fuller and 
Bayrock (1965), and Harington (1984). B. b. bison 
on the other hand is morphometrically more 
remote from B. b. occidentalis and possesses 
more highly derived external characters (e.g. the 
high degree of development of display hair and 
chromatic variegation). 

The results of the present analyses agree with 
the scenario outlined above i.e. 1) the conspeci- 
ficity of Eurasian and North American Holocene 
bison, 2) the dual origin of modern North Ameri- 
can bison from a Beringian element (B. b. occiden- 
talis) and an earlier Nearctic line (B. b. antiquus), 
and 3) subsequent evolution of a grassland- 
adapted form (B. b. bison) and a more conserva- 
tive parkland form (B. b. athabascae), retaining 
more of the ancestral characters. 

Contrary to the interpretation of other North 
American investigators (e.g. Skinner and Kaisen 
1947, Guthrie 1970, Fuller and Bayrock 1965, 


Wilson 1975), McDonald (1981) postulated a 
phylogenesis of modern North American bison 
from an autochthonous Nearctic line (B. antiquus, 
of his classification) which he believes may also 
have given rise to B. b. bonasus. Flerov (1979) 
stands apart from both schools of thought, deriv- 
ing bonasus and athabascae from Palearctic 
priscus, and bison from the Nearctic autochthon 
antiquus. As to the origin of B. b. athabascae, 
McDonald (1981) suggests two alternatives. Body 
size could have been 1) retained from the ancestor 
(his B. antiquus) or 2) resulted from more recent 
adaptive differentiation from B. b. bison. He 
regards the second alternative as probably ‘‘more 
important over time in influencing the nature of 
B. b. athabascae’’ (McDonald 1981, p. 261). As 
we have seen in the preceding analyses more than 
just size is involved. B. b. athabascae differs from 
B. b. bison (and resembles B. b. occidentalis) in 
a number of characters other than size, which we 
may assume to be genetically determined. Con- 
sidering all these characters together, the view that 
B. b. athabascae is a direct and little differentiated 
descendant of B. b. occidentalis is more probably 
correct. 


4. The Nature and Maintenance of Subspecific 
Differences in North American Bison 


As noted in the introduction some regard the 
differences between B. b. athabascae and B. b. 
bison ‘‘as but a superficial divergence springing 
from environment’’ (Soper 1941: p. 375). The idea 
that the observed differences between wood and 
plains bison are a result of environmental induc- 
tion is not new and still finds its adherents. It is 
difficult to separate ecophenotypic variation from 
that under genetic control, without resorting to 
experimental study of the effects of environmental 
variables on morphology. There are, however, a 
few lines of evidence that shed light on the ques- 
tion of whether we are dealing with ecophenotypic 
or genetic variation. 

1. NR bison were transplanted from their origi- 
nal habitat to the Mackenzie Bison Reserve in 1963 
and to Elk Island National Park in 1965, two very 
different environments. Specimens obtained 
14 years later from Elk Island National Park and 
19 years later from the Mackenzie Bison Reserve 
did not differ appreciably from one another, nor 
did they differ appreciably from those specimens 
collected at the original sites in Wood Buffalo 
National Park. Moreover, NR bison from Elk 
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Island National Park do not show to this day any 
evidence of approaching the plains bison from Elk 
Island National Park north of the highway, in the 
morphological characters evaluated, although 
both herds exist under nearly identical environ- 
mental conditions. These observations suggest 
genetic rather than ecophenotypic variation. 

2. The morphometric similarity of NR speci- 
mens, drawn from three widely spaced localities, 
to those of B. b. athabascae and B. b. occiden- 
talis (samples separated in time and space) also 
suggest genetic control rather than environmental 
induction of the observed differences. Likewise, 
the morphometric similarity of AP bison to his- 
toric plains bison, the former living under very 
different conditions than its ancestors, leads to the 
same conclusion. 

3. One can gain some insight into the plasticity 
of the bison phenotype by considering the effects 
of captivity. Plains bison have been kept in zoos 
for many generations. None of the characters 
assessed as characterizing the plains bison has been 
altered significantly — except perhaps for size and 
weight. Captive animals free from stress and on 
a high plane of nutrition, may attain larger size 
and greater weight than their wild counterparts. 
For example, a breeding bull of the Catskill Game 
Farm weighed 1035 kg at age 12 (Heck 1968). That 
weight greatly exceeded the upper limits observed 
in herds living under semi-natural conditions in 
the Wichita Mountains Wildlife Reserve (Halloran 
1960) and in Elk Island National Park. We are 
again led to the conclusion that the characters 
which distinguish wood bison from plains bison 
are genetically controlled, are not ecophenotypic 
variations, and constitute good subspecific 
characters. 

Subspecific distinctness was probably main- 
tained by a combination of ecological and 
behavioral factors as well as by natural selection. 

1. Ecological factors minimized interbreeding 
between the two subspecies, because populations 
were largely spatially isolated by unfavorable 
habitats. The distribution of B. b. athabascae 
within their range was probably highly clumped, 
with demes centered in areas of favorable habitat 
and absent from areas of unsuitable habitat (e.g. 
muskeg). An intervening belt of heavy boreal for- 
est acted as a natural barrier to contact with B. b. 
bison to the south. The barrier, although not 
impenetrable, must have at least minimized con- 
tact. The distribution of B. b. bison populations 
was probably less clumped because of the more 


continuous nature of suitable habitat. 

2. Behavioral factors minimized interbreeding, 
because different traditions tended to separate 
populations during the reproductive season. Bison 
demes tend to occupy relatively restricted home 
ranges of which they possess intimate knowledge 
with regard to resources and their seasonal avail- 
ability. This leads to seasonal movements from 
certain parts of the home range to others. These 
seasonal movements are regular and tend to 
remain stable as long as no major climatic or other 
environmental changes disrupt them. The presence 
of long established and well worn bison trails used 
by generations of bison testify to this (Soper 1941). 
Evidence for such seasonal movements in B. Db. 
bison populations occupying the Canadian prairies 
is presented by Moodie and Ray (1976) and 
Morgan (1980). They documented annual move- 
ments into the boreal parklands in winter and back 
to the prairies in spring and summer. Seasonal 
movements from winter ranges to summer ranges 
and back are also known from Wood Buffalo 
National Park (Soper 1941) and Yellowstone 
National Park (Meagher 1973). Free ranging Euro- 
pean bison, as well, are reported to move 
seasonally (Heptner et al. 1966, Baskin 1979). I 
have no information on the seasonal movements 
of historical populations of B. b. athabascae, but 
we may assume that they did not behave differ- 
ently from free-ranging bison today. Their sum- 
mer ranges were probably remote from those of 
the nearest plains bison. If so, there would have 
been little or no contact between the two subspe- 
cies during the breeding season; an effective 
barrier to gene flow. 

3. Differential directional selection, imposed by 
the diverse habitats of the two subspecies may 
have played a role in maintaining allele frequen- 
cies and phenotypic distinctness. Unlike environ- 
ments can exert a strong selective pressures on 
natural populations. The resulting genetic differ- 
entiation in adjoining populations may confer a 
greater degree of fitness to each, allowing them 
to survive better under the respective specific 
environmental conditions of the different habitats. 
The different selective regimes would thus offset 
any effects of gene flow that may occur. 

The interaction among ecological and 
behavioral factors, gene flow and natural selec- 
tion could thus account for the maintenance of the 
distinctness of B. b. bison and B. b. athabascae. 
Similar mechanisms still operate and can be ana- 
lyzed in other large ungulates, like caribou, (e.g. 
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Rangifer tarandus caribou, R. t. groenlandicus 
and R. ¢. pearyi) which occupy heterogeneous 
environments and have developed very different 
traditions in response to them. 

McDonald (1981) suggests that hunting by 
Indians, which according to him appeared to have 
taken place primarily along the margins of the 
main range, may have been the most influential 
factor separating the athabascae and bison gene 
pools by preventing the peripheral expansion of 
the range of B. b. bison. It is difficult to conceive 
how this could have been achieved given the 
human population densities, which must have 
been quite low in early historic time, and the state 
of contemporary hunting technology. Probably 
human predation in early historic time could not 
have presented a significant barrier to expansion. 


5. Taxonomic Status of Present Remnant 
Populations in Canada. 


Craniometric comparison of AP bison with histor- 
ical populations of B. bison show them to be 
highly similar to historical populations from the 
Canadian prairies. In pelage characters and out- 
line shape AP bison also resemble historical Cana- 
dian prairie bison, as a comparison with photo- 
graphs of an original bison bull (Anonymous 
1906, Lothian 1981) from the Canadian prairies 
showed. The close resemblance of AP bison to 
historical bison from the Canadian prairies is not 
altogether surprising. The origin of the AP herd 
can be traced to stock from the area near the 
Alberta-Montana border and Saskatchewan with 
some admixture of southern stock (Texas, New 
Mexico, Kansas and Nebraska) (Dary 1974, 
Lothian 1981). 

The multivariate morphometric analyses of 
skeletal characters and the comparison of pelage 
characters and outline shape of NR specimens with 
those of B. b. athabascae (insofar as sufficient 
information on homologous characters was avail- 
able) and B. b. bison (AP) specimens show them 
to be somewhat intermediate, but much closer to 
B. b. athabascae. The morphological evidence 
thus suggests that some gene flow occurred 
between the population in the Nyarling-Buffalo 
Lake area and populations elsewhere in the park, 
with B. b. bison genes in their gene pools. Pheno- 
typic similarity cannot be translated into precise 
genetic terms. Morphometric analysis for exam- 
ple measures the effects of many gene loci, but 
how many alleles are responsible for the observed 


divergence is not known. To complicate the pic- 
ture further, non-genetic sources of variation often 
contribute to the results. It is therefore difficult 
to evaluate the degree of hybridization that may 
have occurred. If the differences in the expression 
of frontal display hair and chaps are based on a 
one-locus, two-allele system as suggested earlier, 
the pelage-character analysis suggests a possible 
admixture of AP genes of approximately 5%. 

The conclusion that gene flow occurred is con- 
trary to the original assumption of isolation by 
distance and intervening ‘‘inhospitable and unpro- 
ductive terrain’’ (Banfield and Novakowski 1960). 
Circumstantial evidence from other sources sug- 
gests that the isolation was less than complete. 
Soper (1941) stated that the country between the 
Ninishith Hills and the upper Buffalo River and 
the Caribou Mountains is swampy, except for 
some sandy ridges, inaccessible in summer and 
shunned by bison in winter. However, he mentions 
and illustrates a bison trail which followed the 
Ninishith Hills northwest to the Nyarling River. 
Soper’s informant, Sousi Marie, stated that he saw 
Wainwright bison as far north as the Klewi River 
as early as 1929 and ‘‘that trails, progressively 
dimmer, exist. . . . to the Nyarling.’’ Harington 
(1977) reports an observation by L.P.E. 
Choquette, of trails connecting the Nyarling River 
area with the range of mixed herds, confirming 
the information given by Soper (1941). The pres- 
ence of brucellosis and tuberculosis in the Nyarling 
herd (Novakowski and Stevens 1965) — diseases 
believed to have been introduced into Wood 
Buffalo National Park from Wainwright — also 
point to contact. 

It is unfortunate that too few specimens from 
herds occupying different areas of Wood Buffalo 
National Park were available for this study. A 
morphometric analysis of such material would 
help to interpret the impact of the introductions 
of plains bison on the phenotypic variation of the 
different demes inside and outside the Park. 

If the population estimate of 1500 B. b. atha- 
bascae in 1922 was correct (Graham 1923), the 
greater abundance of B. b. bison in the population 
in the late 1920s (6673) suggests a proportionately 
larger contribution to the genetic make up of the 
present population in Wood Buffalo National 
Park. However, two factors may have influenced 
the degree and extent of hybridization, the demo- 
graphic composition of the population of the two 
taxa and the extent of movement among demes. 
Most of the introduced plains bison were female 
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and young (mostly yearlings and two-year-olds, 
only 332 were three-year-olds). It is therefore 
probable that during the first few years following 
the initial introductions, athabascae bulls were 
primarily involved in breeding. If this was the case, 
the genetic contribution by B. b. athabascae would 
be greater than the gross figures would suggest. 
Release sites of the introduced plains bison were 
all located along the eastern park boundary below 
the confluence of the Peace and Slave rivers. 
Depending on the extent of movement and ex- 
change of individuals among different local popu- 
lations, hybridization may have affected different 
herds (e.g. Lake Claire, Hay Camp and Grand 
Detour herds) to varying extents. If movement 
between demes was limited, the Hay Camp and 
Lake Claire Herd would be expected to show the 
greatest similarity to AP stock. Further work is 
needed to shed light on these problems. 

A comparison of blood characteristics of NR 
bison with those of five different herds of plains 
bison (Peden and Kraay 1979) showed that differ- 
ences among the plains bison herds were at least 
as great as those between NR bison and some of 
the plains herds. The relationship among herds 
based on frequencies of red blood cell antigens are 
not congruent with those based on morphology, 
and did not provide evidence for the subspecific 
distinction. As most of the present herds origi- 
nated from small founder populations, the 
founder effect may have had an important 
influence on the observed frequencies. The lack 
of agreement between morphology and red blood 
cell antigens reflects the fact that different charac- 
ter sets may represent different adaptations 
(assuming that observed differences have adaptive 
significance) subject to different selection pres- 
sures. Blood groups may not have evolved at the 
same rate as morphology, as the high degree of 
similarity of cattle and bison with respect to blood 
groups appears to indicate. Each of the nine blood 
group systems found in cattle has a known homo- 
logue in bison (Stormont et al. 1961). Blood group 
frequencies can be altered relatively easily within 
a species judging from the divergent frequencies 
found in different breeds of cattle. Morphology 
of the body as a whole, on the other hand, appears 
to be linked to the basic ground plan of a species, 
which is probably under complex polygenic con- 
trol and less subject to rapid change. 

Karyotypes of NR and plains bison were found 
to be identical (Ying and Peden 1977), and highly 
similar to that of the European bison (Orlov 


and Chudinovskaya 1979). All possess a com- 
plement of 60 chromosomes. The only difference 
between European and American bison is in the 
Y chromosome, which is submetacentric in the 
former and usually acrocentric in the latter. 
Variation from submetacentric to acrocentric Y 
chromosomes appears to be present in American 
bison (Wurster and Benirschke 1968, Hsu and 
Benirschke 1969, Ying and Peden 1977), but there 
is no information on the frequency of occurrence 
of the two types. The karyotypes of Bison and Bos 
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are highly similar (Orlov and Chudinovskaya 
1979), representing thus another character set of 
low interspecific divergence. 

The results of the present analyses demonstrate 
the close proximity of NR bison to the original 
B. b. athabascae, and warrant their inclusion in 
that taxon, in spite of the apparent slight introgres- 
sion of B. b. bison (AP) genes. The WBNP popu- 
lations could be designated as intergrades (B. b. 
bison X athabascae) until they can be analysed 
in greater detail. 
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Appendix 1 


List of specimens used 


Samples representing remnant populations 

B. b. bison AP males (total 93) — Alberta — Elk Island 
National Park (NMC 6010, 21895, 27504, 27508, 27510, 27511, 
45416, 51387, 51389, 51395; UCYS1, $2, S5, S7, S8, S9, S11, 
Se. O13, S15, SI6, S17, $18, SI9, S21, S22, 824,527, 529, 
S32, S33, S34, S38, S40, S41, S44, S45, S46, S48, S49, S51, 
S52, S55, S56, S59, S61, S63, S64, S66, S68, S70, S72, S74, 
S75, S77, S78, S79, S80, S82, S84, S86, S87, S90, S92, S94, 
$95, S98, S99, S101, S102, S103, S106, S108, S110, S111, S118, 
Seo, 5126, 5127, $128, $130, S131, S132;3133, S134, S135, 
S136, S142, S146, S147). Buffalo National Park, Wainwright 
(NMC 5552). Manitoba — Riding Mountain National Park 
(MMMN 8055, 8554). 

AP females (total 52) — Alberta — Elk Island National Park 
(NMC 6009, 27502, 27503, 27505, 27507, 51388, 51390, 51391, 
51392, 51393, 51394, PMA 66-3-2, 69-35-1, VAMP K980-12, 
UCY S6, S10, S23, S25, S26, S28, S30, S39, S42, S43, S47, 
B51, 553, 557, 967, S69, S76, S81, $83, S85, S88, $93, $96, 
ede, 5110, S114, 5115, 5120, S121, S122, SiZ9, S137, 
$138, S140, S141, S143, S145). 


B. b. bison x B. b. athabascae — WBNP males (total 11) 
— Alberta — Wood Buffalo National Park, Grand Detour 
(CWS Edmonton, unnumbered), Government Hay Camp 
(USNM 263390), Pine Lake (USSC 2330), Lake St. Clair (NMC 
48087, 48089, 48091, 48092, 48093, 48096, 48099, 48100). 

WBNP Females (total 13) — Alberta — Wood Buffalo 
National Park, Sweetgrass (NMC 51399, 51400, 51401, 51428), 
Murdock Creek (NMC 12091), Lake St. Clair (NMC 48084, 
48085, 48086, 48088, 48090, 48095, 48097, 48098). 


B. bison ssp? — SE males (total 6) — Alberta — Wood 
Buffalo National Park (AMNH 86950). Wood Buffalo 
National Park, Salt River (AMNH 130171), Salt Lick, 5 Acre 
Wallow (AMNH 98953), 16 km from Slave River (AMNH 
98228, 98229), Salt Plain (AMNH 98958). 

SE females (total 3) — Alberta — Wood Buffalo National 
Park, Salt Plains (AMNH 98954, 98957, 98959). 


B. b. athabascae — NR males (total 35) — Alberta — Elk 
Island National Park (NMC p. 32628, NMC 45622, 45623, 
45624, 45625, 45626, 45627, 45628, 45629, 45630, 45631, 45632, 
51371, 51374, 51402, 51422, 51423, 51424, PMA 73.82.1, 
74.86.1, 75.23.1. Northwest Territories — Wood Buffalo 
National Park, Needle Lake area (NMC 24026, 24027, 39875, 
51376, 51377, 51378, 51379, 51380), Copp Lake (NMC 51403); 
Mackenzie Bison Sanctuary (NMC 51408, 51409, 51412, 51413, 
51417). 

NR females (total 19) — A/berta — Elk Island National Park 
(NMC 51372, 51405, 51406, 51407, 51425). Northwest Terri- 
tories — Wood Buffalo National Park, Frying Pan Lake (VMC 
24028), Needle Lake area (NMC 39876, 51381, 51382, 51383, 
51384, 51385, 51386, 51427); Mackenzie Bison Sanctuary 
(NMC 51410, 51411, 51414, 51415, 51416). 
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B. b. bonasus — males (total 23) — BM(NH) 45.10.13.2, 
1946.12.17.1; BZM 8030, 20466, 22072, 22077, 33832, 33833, 
33824, 33863, 41709, 44213, 46548, 91474, ZMA 1400, 1404, 
5290.16.037, 18.404; Wien 1311, 1558, 2765, 5536, 7056. 

Females (total 11) — BM(NA) 92.3.15.2; BZM 8032, 14798, 
22078, 33837, 33838, 45287; IH 4623; MMMN 8596, 1937.153; 
ZMA 1335. 


B. b. caucasicus — males (total 4) — BZM 33861; Wien 1594, 
1596, 2320. 
Females (total 2) — BZM 41673; Wien 2322. 


Samples representing historical populations 
B. b. athabascae — males (total 19) — Alaska — valley of 
Tanana River (USNM 223292). Alberta: Wood Buffalo 
National Park (AMNH 73615, PMA 68-31-1, 68-31-2, 68-31- 
4, UAMZ p626); near Fort Smith (NMC 4538) west of Fort 
Smith (USNM 177632); Lying Wood Mountain (NMC 625); 
80.5 km SW of Smith Landing (USNM 172689); between 
Athabasca and Great Slave Lake (MVCZ 24017); Wood Buffalo 
National Park, Salt Prairie (WMC 8755), Murdock Creek, 
(NMC 11436), Blind Bull Lake (NMC 11435). Northwest 
Territories: Great Slave Lake (AMNH 18163, MCZ 21933); 
80.5 km SW Fort Resolution (NMC 299); near Fort Providence 
(USNM 177630); Yukon: Pelly River near Selkirk (USNM 
5515). 

Females (total 1) — Alberta: Pine Lake, Wood Buffalo 
National Park (NMC 10405). 


B. b. bison — CP males (total 58) — Alberta — Drumheller 
(NMC 1803, PMA 80-10-1); junction of Bow and Belly Rivers 
(NMC 5591); Little Sandhill Creek near Red Deer River (VMC 
5592); Beaver Crossing (NMC 51373); Trout Creek (UCY A: 
05-825); Big Hill Springs (UCY A(G); Jenkin Buffalo Jump 
(UCY A (I)); S. of Spruce Grove (Alta. Beach Mus. unnum- 
bered) Alexander Indian Res. (Alta. Beach Mus. unnumbered); 
Alberta Beach (Alta. Beach Mus. unnumbered); Fort 
Assiniboine (Barrhead Mus. unnumbered); Pembina Flats 
(Sangudo Mus. 3 unnumbered); W. side Cooking Lake (Tofield 
Mus. unnumbered); 9.7 km W. Tofield (Tofield Mus. 2 
unnumbered); Beaverhill Lake (Tofield Mus. unnumbered); 
Buffalo Lake, 7 km SE Bashaw (priv. coll.); Medicine Hat 
(priv. coll.); Sundre (priv. coll.); Fort Saskatchewan (priv. 
coll.); South Ghost River (priv. coll.); Bonnyville (priv. coll.); 
Hairy Hill (priv. coll.). Manitoba: near Lake Winnipegosis 
(NMC 9564); Gypsumville (MMMN 8073); McArthur Lake, 
Riding Mountain National Park (NMC 33386). Saskatchewan: 
(NMC_ 10345); bank of South Saskatchewan River near 
Saskatoon (US 2299); near Chaplin (NMC 9565); S. of 
Battleford (NMC 9566, 9567); near Prince Albert (NMC 8833); 
Old Wives Lake (NMC 51426, UCY AC(J)); shore Last 
Mountain Lake near Silton (SMNH p632); NW of Regina 
(SMNH p537); Killdeer (SMNH 10922); between Bolney and 
Maidstone (PMA 79-26-1); Osler (17 priv. coll. unnumbered). 


CRM males (total 53) — Alberta: Between Red Deer and 
Panther River S of Yaha Tinda Ranch (NMC 33395); Yaha 
Tinda Ranch (NMC 33372, 33363); Buffalo Lakes (VMC 
33355); Bowden (NMC 33361); Smoky River (NMC 33393); 
Banff National Park, Stoney Creek (NMC 33365); Red Deer 
River, 12.9 km W Scotch Camp (NMC 33354, 33384); Baker 
Creek (NMC 33378), Panther River (NMC 33385), Tyrrell 
Creek Flats (VMC 33399), Tyrrell Creek (NMC 33371), Indian 
Head Warden Station (NMC 33389), Glacier Lake (VMC 
33379), Bow River, 8 km W of Banff (NMC 33391), 9.7 km 
S of Dormer Cabin (NMC 33370), 19.3 km up Spray River 
(NMC 51396); Jasper National Park, Sulphur Springs (VMC 
33358), Brazeau District (NMC 33356, 33374, 33390), Brazeau 
Falls (NMC 33357), Jasper Park Lodge (NMC 33406), Moose- 
horn River (NMC 33376), Maligne Lake (NMC 33377), Isaac 
Creek (NMC 33362), lower end Maligne Range (NMC 33381), 
Willow Creek (NMC 33368, 33375); Waterton Lakes National 
Park (NMC 33360, 33366, 33401, 33405, 36173, 51397, 51398); 
Waterton Lakes National Park, Twin Butte (priv. coll. 3 
unnumbered), Spring Creek Hatchery (NMC 33340, 33392), 
N Fork Trail Creek (NMC 33364), Canyon Creek Church 
Camp (NMC 51370), Pass Creek Flats (NMC 33403), north 
side Trail Creek (NMC 33359), Horseshoe Basin (NMC 33398), 
mouth of Pass Creek (NMC 33397, 33380), Waterton River 
(NMC 33404, 33396), Indian Creek Badlands (NMC 33382, 
33402). 

CRM females (total 4) — Alberta: Mount Sarbach, Banff 
National Park (NMC 33388); Jasper National Park (NMC 
33373); Jasper National Park (NMC 33367), Willow Creek 
(NMC 33369). 

PR males (total 24) — Alberta: Kinuso (NMC 51375); Man- 
ning (NMC 51418, 51419); 16 km N DeBolt (DeBolt Mus. 
unnumbered); Brownvale (Brownvale Mus. unnumbered). 
Eureka River, 48.3 km NW Fairview (Fairview Mus. unnum- 
bered); N of Enilda (High Prairie Mus. X38); W of Triangle 
(High Prairie Mus. X40-1-2); 38.6 km S of High Prairie (High 
Prairie Mus. X56), High Prairie area (High Prairie Mus. X63); 
1.6 km N of Enilda (High Prairie Mus. 973-101); 9.7 km S of 
town of Peace River (Pioneer Mus. 73-557); 80.5 km from town 
of Peace River (Pioneer Mus. Grand Prairie 20); Wabasca (priv. 
coll.) Spirit River (priv. coll.); 24.1 km WNW Brownvale (priv. 
coll.), Calling Lake (pers. comm. N. Novakowski). British 
Columbia: Fort St. John area (BCPM 2794, priv. coll. unnum- 
bered); Upper Cache Creek (NMC 51420, priv. coll. unnum- 
bered); Pouce Coupé River, 16 km NE Dawson Creek (Dawson 
Creek Mus. unnumbered); Hudson Hope area (Hudson Hope 
Mus. unnumbered); Rose Prairie (priv. coll.). 

PR females (total 1) — Alberta Torrens Towers (PMA 
69-93-1). 

CGP males (total 45) — Montana: AMNH 5475, 5476, 5480, 
130208, 130209, 16300, 16301, 16303, 16305, 16306, 16307, 
16308, 16311, 16313, 16314, 16315, 16317, 16318, 16321, MCZ 
8377, 34083, USNM 17632, 22664, 22665, 22667, 122667, 
122669, 122670, 122673, 122674, 122676, 122677, 122688, 
122689, 122693, 122695, 122699, 233520, 243298, 274791, 
300324. Wyoming: MCZ 1215, 1216, 1218, 7895. 
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CGP females (total 22) — Kansas: MCZ 105. Montana: 
AMNH 130210, 130211, 130212, 16322, 16323, USNM 2841, 
2842, 38302, 38306, 122669, 122672, 122685, 179087, 300324. 
South Dakota: AMNH SD-2180, SD-2181, SD-2182, SD-2183, 
SD-2184. Wyoming: MCZ 1217, UWYO 3081. 

USRM males (total 16) — Colorado: DMNHC 1846, USGS 
A. Montana: UM G0295, USNM 120507, 120512, 221089. 
Oregon: USNM 246529, 249842, 249844, 249846, 249894, 
249895, 250092. Utah: USNM 1442. Wyoming: AMNH 
139937, USNM 168816. 

USRM females (total 13) — Oregon: USNM 249841, 249843, 
249845, 249847, 249848, 249849, 249850, 250089, 250090, 
250091, 250093,.250094, 250095. 

SGP males (total 15) — New Mexico: USNM P(A). Okla- 
homa: TMM A. Texas: PPHM B, C, D, F, 630a, 76 2-1, 915- 
a, 1173-1, 1482-1, TMM 139, TTU A13212, USNM 126620, 
126621. 

EUS males (total 2) — Kentucky: MCZ 2047, 2050. 


Samples representing earlier Holocene and late Pleistocene 
populations 
B. b. occidentalis — males (total 43) — Alaska: USNM p4157. 
Alberta: NMC 12442, 17530; PMA p81-42-1, p68-2-150, p68- 
2-164, p68-2-195, p-68-2-624, p-69-17-19, p74-16-2; Fort 
Saskatchewan priv. coll. 2 unnumbered; near Calgary priv. coll. 
unnumbered; Banff Archives, unnumbered, pers. comm. C.R. 
Harington (1984); Fort Saskatchewan priv. coll. unnumbered; 
Pioneer Mus. Grande Prairie 18, 1; Spirit River NMC 51421; 
Iowa: ISUMZ 38959; SMP 21-67-2, 43-51-2, 155-56-2, 156-56-2 
(A). Manitoba: NMC 2813. Minnesota: MMNH 198, 450, 723, 
738, USNM p10541, p10545. Nebraska: FMNH p26249. North 
Dakota: NDSU G166. Saskatchewan: SMNH p1240a, p1240c. 
South Dakota: ISUMZ p199. Texas: DMNHT 629, TMM 725, 
TTU A1975/18. Wisconsin: MMNH 603A, 603E. Yukon: 
MNC p2242, 1733, USNM pS5513. 

Females (total 6) — Alberta: PMA p68-2-148. Minnesota: 
SMM p67-1-34. Wyoming: UWYO NC2560, p69-20-1. Califor- 
nia: LACM RLB913, RLB1007. 


B. b. antiquus — males (total 4) — Alberta: NMC p12442, 
PMA p69-20-1 — California: LACM RLB913, RLB1007. 

Females (total 6) — California: CM 10195, LACM HC1245, 
1246, 1255, YPM 10602 — Minnesota: USNM p10544. 


B. b. priscus — males (total 3) — Alaska: AMNH A685- 
3249, AMNH AFAM30623, AMNH AFAM46885. 


Other: Eastern Siberia — Moscow Pal. Mus. 835-624-39 
(male), Netherlands — ZMA 4388 (male). 


Appendix 2 Discrimination of Crania of Adult Male B. b. bison and B. b. athabascae 
Discriminant Function Analyses Using Different Numbers of Variables 


Using nine characters 


Discriminant athabascae (N = 12) bison (N=75) Score contribution 
Character Weights character means character means to distance (D2) 
VD 0.391508 94.06 74.47 7.668 
TD 0.090249 98.72 79.15 1.766 
WHCO 0.201389 295.67 265.35 6.105 
OC 0.179051 186.17 F51..51 6.205 
CW -0.110337 127.38 123.66 -0.410 
CLEC 0.073058 292.67 241.96 3.705 
TB -0.093566 Ad ea 179.07 -3.191 
GPW -0.051882 358.00 320.03 -1.978 
SHTT 0.000866 700.75 604.57 0.083 
Dividing point between groups 98 .046* 
athabascae 
centroid 108.027 
range 98.270 - 114.749 
bison 
centroid 88.066 
range 78.448 - 98.219 
D2 athabascae — bison 19.961 
Standardized normal deviate 2233 


Probability of incorrectly assigning an unknown specimen = 0.01 


* The dividing pounit is shifted from the mid point toward the less numerous group by an amount equal to (log R)/D, R is 
the ratio of the more numerous to the less numerous group (Davies 1971). 


Using eight characters 


Discriminant athabascae (N = 12) bison (N=74) Score contribution 

Character Weights character means character means to distance (D7) 
VD 0.485495 94.06 74.47 9.443 
ID 0.010876 98.72 79.39 0.210 
WHCO 0.099340 295.67 266.22 2.926 
OC 0.117815 186.17 153.28 3.874 
CW -0.109679 127.38 123.82 -0.389 
GELC 0.092826 292.67 241.20 4.777 
TB -0.107834 213.47 180.36 -3.537 
GPW 0.050245 358.00 321.50 1.834 
Dividing point between groups 96.767 
athabascae centroid 106.241 

range 98.522 - 113.075 
bison centroid 87.104 

range 76.112 - 99.227 
D2 athabascae — bison 19.137 
Standardized normal deviate 2.181 


Probability of incorrectly assigning an unknown specimen = 0.01 
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Using seven characters 


Discriminant athabascae (N = 12) bison (N=75) Score contribution 

Character Weights character means character means to distance (D2) 
VD 0.510904 94.06 74.68 9.900 
TD 0.018877 98.72 79.35 0.366 
WHCO 0.137628 295.67 266.33 4.037 
OC 0.129480 186.17 153.13 4.277 
CW -0.119818 127.38 123.81 -0.427 
CELE 0.090518 292.67 240.99 4.678 
TB -0.111933 21317 180.45 -3.662 
Dividing point between groups 92.595 
athabascae centroid 102.084 

range 93.454 - 108.753 
bison centroid 82.914 

range 73.065 - 95.103 
D2 athabascae — bison 19.170 
Standardized normal deviate 2.189 


Probability of incorrectly assigning an unknown specimen = 0.01 


Using five characters 


Discriminant athabascae (N = 12) bison (N = 88) Score contribution 

Character Weights character means character means to distance (D2) 
VD 0.318622 94.06 74.48 6.239 
TD 0.058985 98.72 79.30 1.146 
WHCO 0.107565 295.67 266.06 3.185 
OC 0.142588 186.17 152.17 4.847 
CW -0.071918 127.38 123:92 -0.263 
Dividing point between groups 77.534 
athabascae centroid 84.979 

range 76.810 - 91.577 
bison centroid 69.826 

range 61.312 - 83.542 
D2 athabascae — bison 15.154 
Standardized normal deviate 1.946 


Probability of incorrect identification of an unknown specimen = 0.02 
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Using four characters 


Discriminant athabascae (N = 12) bison (N=75) Score contribution 

Character Weights character means character means to distance (D2) 
VD 0.407611 94.06 7.19 7.714 
TD -0.014046 98.72 79.79 -0.266 
WHCO 0.064110 295.67 267.64 1.797 
oC 0.099538 186.17 154.05 3.196 
Dividing point between groups 68.365 
athabascae centroid 74.439 

range 66.506 - 79.884 
bison centroid 61.997 

range 53.495 - 71.960 
D2 athabascae — bison 12.441 
Standardized normal deviate 1.764 


Probability of incorrect identification of an unknown specimen =~ 0.04 


Using three characters 


Discriminant athabascae (N = 12) bison (N = 97) Score contribution 

Character Weights character means character means to distance (D2) 
VD 0.464413 94.06 74.62 9.028 
TD -0.050373 98.72 79.43 -0.971 
WHCO 0.126151 295.67 266.74 3.649 
Dividing point between groups 70.275 
athabascae centroid 76.000 

range 68.637 - 81.091 
bison centroid 64.302 

range 57.825 - 73.413 
D2 athabascae — bison 11.698 
Standardized normal deviate 1.710 


Probability of incorrect identification of an unknown specimen = 0.04 
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